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\Q ' ABSTRACT 

<N 

Context. Searching for extrasolar planets through radial velocity measurements relies on the stability of stellar photospheres. Several 
phenomena are known to affect line profiles in solar-type stars, among which stellar oscillations, granulation and magnetic activity 
(**V| through spots, plages and activity cycles. 

Xf\ , Aims. We aim at characterizing the statistical properties of magnetic activity cycles, and studying their impact on spectroscopic 

• . measurements such as radial velocities, line bisectors and line shapes. 

i-^ ■ Methods. We use data from the HARPS high-precision planet-search sample comprising 304 FGK stars followed over about 7 

»-^H' years. We obtain high-precision Ca II H&K chromospheric activity measurements and convert them to R' HK indices using an updated 

*-*. . calibration taking into account stellar metallicity. We study R' HK variability as a function of time and search for possible correlations 

5^ , with radial velocities and line shape parameters. 

Results. The obtained long-term precision of ~0.35% on S-index measurements is about 3 times better than the canonical Mt Wilson 

survey, which opens new possibilities to characterize stellar activity. We classify stars according to the magnitude and timescale of 

the Ca II H&K variability, and identify activity cycles whenever possible. We find that 39±8% of old solar-type stars in the solar 

neighborhood do not show any activity cycles (or only very weak ones), while 61+8% do have one. Non-cycling stars are almost only 

found among G dwarfs and at mean activity levels log^^ < -4.95. Magnetic cycle amplitude generally decreases with decreasing 

^- ' activity level. A significant fraction of stars exhibit small variations in radial velocities and line shape parameters that are correlated 

l/") | with activity cycles. The sensitivity of radial velocities to magnetic cycles increases towards hotter stars, while late K dwarfs are 

C*\| . almost insensitive. 

Conclusions. Activity cycles do induce long-period, low-amplitude radial velocity variations, at levels up to ~25 ms _1 . Caution 
l/~) ' is therefore mandatory when searching for long-period exoplanets. However, these effects can be corrected to high precision by 

1^* detrending the radial velocity data using simultaneous measurements of Ca II H&K flux and line shape parameters. 

Key words. Planetary systems - Stars: activity - Stars: chromospheres - Line: profiles - Techniques: radial velocities - Techniques: 
spectroscopic 



- * — i 

X 



1. Introduction radius to properly derive exoplanet parameters. The behavior 

of stellar photospheres also plays a crucial ro le: since the dis- 

b , The advent of large-scale Doppler surveys to search for extra- covery of 51 p eg b dMavor & Ouelozl fl995h . questions have 

solar planets around FGKM stars in the solar neighborhood has adsen about the eflfects of steUar photospheric "jitter" on the 

produced an impressive body of high-resolution spectroscopic detection and characterization of exoplanets (e.g. lBaliunas et"aTI 

data extending over the past 20 years or so. Precise radial ve- [T997K The outer conve ctive envelope that is present in solar- 

locities are obviously the main products that are derived from type stars gives rise t0 several phenomena that can have an im- 

these data, but more generally the existing libraries of spectra pact on derived ex0 planet properties, through variations in mea- 

make it possible to study in details many properties of solar- sured d i sk . aver aged radial velocity and photometric flux. In or- 

type stars. The knowledge of extrasolar planets has always been der of increas i ng ti mescales, one can mention p-mode oscilla- 

intimately related to the knowledge of their parent stars. An tions (Christensen-Dalsgaard 2004), granulation and supergran- 

important example is the need for precise fundamental stellar ulation (H arvevlll985h. m^netic a ctivity ind ucing surface inho- 

parameters such as effective temperature, metallicity, mass and mogen eities rotating with the star (ISaar & Donahuelfl997l) . and 

finally magnetic activity cycl es over timescales of years and 



Based on observations made with the HARPS instrument on the decades (BaUunas et al. 1995]). From a stellar physics point of 



ESO 3.6-m telescope at La Silla Observatory (Chile), under program v j eWi a q tne se phenomena are obviously interesting in their own 
IDs 072.C-0488, 183.C-0972, 083.C-1001 and following periods. 
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right since they can reveal the details of stellar internal structure 
and dynamo mechanisms at play in solar-type stars. Moreover, 
the statistical study of stellar activity in a large number solar- 
type stars can also inform us about the past and future evolution 
of our Sun, and its potential impact on Earth climate. 

With the ever-increasing precision of radial velocity (RV) 
and transit observations of exoplanets, the need to better un- 
derstand the host stars has become more acute in recent years. 
In particular, radial veloc ity surveys are now probing the sub- 
ms -1 regime (Pepa l201ll) . a level at which stellar activity per- 
turbations become non-negligible even for quiet stars and must 
be monitored. The study of the effect s of star spots on radia l 



velocities has a l o ng history, see e.g. j Saar & Donahue ( 
Saar et alj (f l998); S antos et all (J2000); Saar & Fischer 



1997); 



(2000) 



Oueloz et al.1 yOOl ). More recently, the field has received closer 



attention again following the discovery of transit ing planets 
around relatively active stars like CoR oT-7 (lOueloz et al.1 12009: 



Bois se et al.ll201 it lHatzes et al.ll201 ll) or M dwarfs like GJ 436 



(Knu tson et al.M2.Ql II) . At the same time, efforts to characterize 
RV jitter in quiet solar-type stars, in particular the Sun itself, 
have also intensified. For example, Dumusque et al. (201 lcb) 
simulate the impact of p-mode oscillations, granulation, super- 
granulation and active regions on RV data, b ased on astero- 
seismic obs ervations of a few so lar-type stars. Lagrange et al.l 
(2010) and M eunier et al.l d2010h use solar irradiance and Ca 
II archival data to infer the level of RV jitter caused by spots 
and plages. While these studies offer important insights into the 
causes of RV jitter, they rely at least partly on unverified as- 
sumptions to relate a given physical phenomenon at the stellar 
surface to the corresponding RV perturbation, as actually mea- 
sured by RV planet-search techniques. The ultimate "truth" in 
this domain must necessarily come from actual measurements 
of precise RVs, photometry and/or spectroscopic diagnostics in 
the form of densely-sampled time series covering the relevant 
timescales. The present paper is an attempt is this direction, fo- 
cusing on the longer timescales. 

Over years and decades, the main cause of variability in 
the Sun is the occurrence of a magnetic cycle with a quasi- 
periodicity of 1 1 years (or 22 years depending on definition). 
The study of this cycle has a long history, starting with its dis- 
covery by S. H. Schwabe in 1843. The question whether a sim- 
ilar phenomenon also occurs in other stars was addressed in the 
1960s by O. Wilson, who started what would become the famous 
Mt Wilson program for measuring chromospheric emission in 
stars of the solar neighborhood (Wilson 19M Il978l) . The Mt 
Wilson spectrophotometers regularly monitored the flux in the 
core of the Ca II H & K lines in the sample stars for more than 
30 years, yielding what is still today the reference database for 
magnetic activity in solar-type stars. Long-term re sults from this 
program were publ i shed in several pape rs (e.g. I Wilson! 119781: 
Dunc an et al.| [l991: Baliuna s et al.lll995l) . which demonstrated 
that activity cycles similar to the solar one are indeed common 
in other stars. 

In the meantime, other stellar activity surveys, and the de- 
velopment of large-scale Doppler exoplanet programs, have also 
produced a large quantity of Ca II H & K observations, and sev- 
eral studies have been conducted to investigate magnetic activ- 
ity levels in solar-type sta r s (e .g. Henr y et alii 199 6: Santos et al. 



2000; Tinnev et al. 



pe stai 
120021 



Wrig ht etal.ll2004t lHall et al.l 120071 



Isaacson & Fischer 2010; Santos et al. 2010). All these "second- 
generation" projects tried to reproduce as closely as possible the 
method used by the Mt Wilson survey to measure Ca II H & K 
emission, i.e. by computing the so-called S index which normal- 
izes the measured Ca II core flux by the flux in two "continuum" 



bandpasses on the blue and red sides of the Ca II lines. The Mt 
Wilson scale, as defined by the S index, has therefore become 
the standard way of measuring chromospheric activity in stars. 

Most projects have actually focused on the mean activity 
level of stars, and not on the time variability of t he ac tivity. 
To the best of our knowledge, only iBaliunas et al.l d!995l) . us- 
ing the Mt Wilson data, have actually performed a systematic 
search for magnetic cycles and derived their main parameters 
such as cycle period. In the present paper we attempt to do a 
similar study to better characterize overall cycle properties in 
older solar-type stars, and compare them to the Mt Wilson sur- 
vey. Then we search for the effects of magnetic cycles on the 
shapes of photospheric lines, with the ultimate goal of disentan- 
gling stellar RV "jitter" from real barycentric motions of the star. 

The idea that magnetic cycles may influence spectral lines 
dates back at least to Dravins ( 1985). The rationale behind this is 
that the convective patterns at the surface of solar-type stars may 
change along the magnetic cycle under the influence of changing 
magnetic field strength. Indeed, the convection is greatly reduced 
in active regions, which causes a decrease in the convective 
blueshift usually exhibited by photospheric lines (see Dravins 
1982; Livingston 1982; Brandt & Solanki 1990; Gray et al. 



1992; Lindegr en & Drav ins 2003; Me unier et al . 2010, and ref- 
erences therein for more detailed studies of these effects). 

A few attempts to measure such an effect in the Sun yielded 
somew hat contradictory results: indeed, Deming & Plymate 
(1 19941) find a peak -to-peak RV va r iation of 28 ms _1 over the 
solar cycle, while [McMillan et al.l (1 19931) obtain constant RVs 
within ~4 ms _I . We note here that the two studies used very 
different spectral lines for their measurements (CO lines in 
the IR vs. atomic UV lines). Measurements of similar RV ef- 
fects in other stars have remained inconclusive to da te (e.g . 
ICampbell et al.lll988b ISaar & Fischerll2000t ISantosetalJ fcolO). 
although Sant os et al.1 Q2010) find clear correlations in several 
stars between magnetic cycles and line shape parameters like bi- 
sector and line width. 

In the present pape r we use high-pr ecision data from the 
HARPS spectrograph ( Mayor et al.ll2003l) . obtained in the con- 
text of a planet-search survey around FGK dwarfs in the solar 
neighborhood. This survey focuses on the search for low-mass 
planets, i.e. Neptunes and super-Earths, which induce RV varia- 
tions of only a few m s" 1 on their parent star (see e.g. Lovis et all 



2006: lMavoretal.l2009l:lLovisetal.l2011l;lPepel2011l) . Potential 
effects of magnetic cycles at this level could induce spurious 
long-term drifts in the data that may perturb the detection of low- 
mass objects at shorter periods due to limited sampling, or even 
mimic the signal of giant planets orbiting at several AUs from the 
star. It is therefore crucial to understand the effects magnetic cy- 
cles may have on our sample stars. I n an ac companying paper to 
the present study, Dumusque et al.l (1201 lal) show several exam- 
ples of planetary systems around stars showing magnetic cycles, 
and discuss how to correct the perturbing effects of these. 

2. Observations and data reduction 

2. 1 . The stellar sample 

The original HARPS FGK high-precision sample comprises 
about 400 FGK stars chosen fro m the volume-lim ited (50 pc) 
CORALIE sample of -1600 stars (lUdrv et al.ll2000T) . The selec- 
tion has been made based on criteria to minimize stellar radial 
velocity jitter and thus maximize the detectability of low-mass 
exoplanets. The main selection criteria were: 1) vsin/ < 3-4 
km s , and 2) exclude all kinds of binary or multiple stars (based 
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Fig. 1. H-R diagram of the 304 stars in the sample. Subgiants are 
shown as open dots. 



on CORALIE measurements). The main effect of this is to gen- 
erate a stellar sample that is biased against young stars, due to the 
v sin i cut. Otherwise, we expect that this sample is still represen- 
tative of old (single) solar-type stars in the solar neighborhood. 

The sample stars have been regularly monitored with 
HARPS since 2003. Single observations usually reach SNR > 
100 per pixel at 550 nm, in order to achieve a radial velocity 
precision of 1 ms _1 or better. Observing cadence varies signif- 
icantly from star to star, depending on radial velocity variabil- 
ity, the presence of planetary systems, and activity level. Active 
stars with log R' HK > -4.7 that could not be eliminated based on 
CORALIE measurements were screened early by HARPS and 
only rarely measured. As of 2011, the number of data points 
per star varies between 1 and more than 200, which makes any 
global analysis of the sample challenging because of the very 
irregular sampling. 

For the purpose of this paper, we further selected stars with at 
least 6 measurements as of May 201 1, and with a total time span 
of observations of at least 1000 days. We finally obtain a sam- 
ple of 304 stars, still representative of old FGK stars in the solar 
neighborhoo d. Fundame n tal pro perties of these stars have been 
obtained by ISousa et all J2008I) using the same HARPS spec- 
tra that we use here. Their analysis yields effective temperature, 
metallicity and surface gravity, while the Hipparcos catalogue 
provides absolute magnitude and luminosity through a bolomet- 
ric correction. Fig. Q] shows a H-R diagram of the 304 stars that 
we selected for this study. They span a range in effective tem- 
perature between 4600 and 6200 K, and a range in metallicity 
[Fe/H] between -0.84 and +0.39. 

We immediately see from Fig. Q] that some stars have al- 
ready evolved off the main sequence, i.e. have become subgiants. 
Because of their lower gravities and higher luminosities, chro- 
mospheric activity measurements in these stars cannot be cali- 
brated in the same way as for main-sequence stars, and we will 
therefore consider them separately in the remainder of this pa- 
per. To separate subgiants from main-sequence stars, we fitted 



the main sequence in Fig. Q] and place a limit at +0.25 dex in 
luminosity (0.63 mag) above the main sequence. According to 
this criterion, there are 20 subgiants among the 304 sample stars, 
shown as open dots in the figure. 

2.2. HARPS CallH&K index 

We compute the HARPS S index by following the original 
Mt Wilson procedu re as closely as possible (see in particular 
iDuncan et alj[l9 9l!). We work on extracted 2D spectra in (order, 
pixel) space, as produced by the standard HARPS pipeline. The 
spectra are wavelength-calibrated, flat-fielded and de-blazed us- 
ing calibration frames acquired at the beginning of each night. 
Their wavelength scale is then transformed to the rest frame 
of the star using the measured stellar radial velocity and the 
barycentric correction. Besides this standard treatment, there are 
two other key steps in the data reduction that are particularly im- 
portant for S index computation and other spectrophotometric 
measurements: 

- Background subtraction: there is a small but non-negligible 
amount of diffuse light in the spectrograph that must be sub- 
tracted to preserve spectrophotometric accuracy and preci- 
sion. In HARPS the diffuse light has two origins: stellar light 
(fiber A) and ThAr light (fiber B) when observing in simul- 
taneous reference mode. Scattered stellar light does not af- 
fect precision to first order since its level is directly pro- 
portional to the total stellar flux entering the spectrograph. 
More problematic is the diffuse light from ThAr because 
it has a constant flux level, independent of stellar flux. As 
a consequence, a clear dependence of the S index on SNR 
is produced, particularly towards low SNR. Correcting the 
background is therefore essential for spectrophotometric ac- 
curacy and precision. We measure the background locally on 
the detector using the inter-order pixels and taking the mode 
of the pixel value distribution as the background value. Then 
the local values are interpolated with splines in the two di- 
rections to generate a global background frame that is sub- 
tracted from the science frame. 

- Contamination by saturated ThAr lines: besides producing 
diffuse light, the ThAr spectrum also contains a number of 
saturated lines that contaminate neighboring stellar spectral 
orders on the detector. To correct this, we developed a new 
spectrum extraction algorit hm based on the usual Home op- 
timal extraction algorithm dHornd[l986l) . We modified the 
algorithm to simultaneously fit the order profile in cross- 
dispersion direction and the contaminating flux coming from 
the neighboring fiber. We obtained a reference ThAr spec- 
trum with fiber B only (no light on fiber A), and use it as a 
template for contamination on fiber A. Only one extra free 
parameter is required in the extraction procedure (a global 
scaling of the contamination level). 

We verified that these improvements removed most of the 
strong systematics as a function of SNR that were affecting Ca 
II H & K measurements before any background subtraction was 
applied. S-index measurements on standard stars show that we 
can now still get reliable results at SNR values of only a few per 
pixel in the continuum at 4000 A. Residual instrumental system- 
atics as a function of SNR are still visible on some stars, and 
their potential causes should be investigated further (e.g. CCD 
non-linearity). They are however negligible for the purpose of 
this paper. 

To compute the S index, we define exactly the same spec- 
tral passbands as for the Mt Wilson HKP-2 spectrophotometer 
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(lDuncanetal]ll99lh . For the Ca II H & K line cores, the bands 
are centered at 3933.664 A (K) and 3968.470 A (H), and have 
a triangular shape with a FWHM of 1 .09 A. For the continuum 
passbands, the V band has a width of 20 A and is centered at 
3901.070 A, while the R band has the same width and is cen- 
tered at 4001.070 A. The Mt Wilson S index is then defined as: 



S — a ■ 



H + K 
R+V' 



(1) 



In this equation, H, K, R and V represent the total flux 
in each passband, while a is a calibration constant. The 
value of a is repor ted as being either fixed at 2.40 or 2.30 
(Duncan et al. 1991), or regularly adjusted following calibration 
dBaliunas et al. 1995). We note here two important facts: 

- The Mt Wilson spectrophotometer is designed in such a way 
that a rapidly rotating slit mask sequentially observes the 
four different channels at high frequency, and exposes the 
H and K channels 8 times longer than the continuum pass- 
bands. 

- When working with 'normal' spectrographs like HARPS, we 
will then get 8 times less flux in the H and K passbands than 
the Mt Wilson spectrophotometer. To be on the Mt Wilson 
S-index scale, we have to multiply the H and K fluxes by a 
and then by 8. 

Now, instead of working with integrated fluxes in each pass- 
band, we could also consider working with mean fluxes per 
wavelength interval, i.e. X - X/AAx, where X represents any 
bandpass and AAx is the effective bandpass width. The Mt 
Wilson S index can then be written as: 



_ AA HK H + K 

a ■ 8 • • — — = a ■ 



1.09 A H + K 



A/1 



rv R+V 



20 A R + V 



H + K 
R + V 



(2) 



We see that, accidentally or not, S is actually equal (or very 
close) to the ratio of the mean fluxes per wavelength interval. 
Indeed, with a - 2.3, the numerical factor is just 1.0028. The 
simplicity of this expression led us to actually use it, i.e. work 
with mean fluxes per wavelength interval instead of integrated 
fluxes. One advantage of this approach is that it helps minimize 
potential edge effects at band boundaries due to finite pixel size 
and errors in wavelength calibration. 

Together with the S value, we also compute the uncertainty 
on S due to photon shot noise through error propagation. 

2.3. Calibration to Mt Wilson scale 

As explained above, we expect that the HARPS S index will al- 
ready be very close to the Mt Wilson scale, so that in principle 
no further calibration should be needed. However, at least two 
aspects make it still necessary: potential scattered light in the Mt 
Wilson spectrophotometer and/or HARPS, and the exact value 
of the calibration constant a. In any case, a simple linear cali- 
bration between the HARPS and Mt Wilson systems should be 
sufficient. 

The problem of calibrating the S scale is that a number 
of 'standard' stars observed with both instruments are needed. 
However, many stars show intrinsic variability because of rota- 
tional modulations and/or magnetic cycles that will induce noise 



Table 1. Refer ence stars for S-ind ex calibration. Mt Wilson val- 
ues are from lBaliunas et aD(ll995h . 



Star 



HD 10700 


0.171 


0.142 


HD 16160 


0.226 


0.194 


HD 23249 


0.137 


0.108 


HD 152391 


0.393 


0.341 


HD 160346 


0.300 


0.250 


HD 216385 


0.142 


0.119 


HD 219834 


0.155 


0.121 



in the calibration procedure. We choose to limit ourselves to a 
small sample of stars that have been well observed at Mt Wilson 
and by HARPS, and whose intrinsic variability remains moder- 
ate. Several stars come from a dedicated HARPS pro gram aim- 
ing at following Mt Wilson stars on the long-term (iSantos et al.l 
2010). Table Q] gives the list of calibrators with their mean Mt 
Wilson and H ARPS values. Mt Wilson values are taken from 
Baliu nas et al.l d 19951) . The calibration stars span a range in S 
values that almost covers the range of values encountered in our 
sample. A linear fit to Mt Wilson and HARPS values yields the 
desired calibration: 



= 1.111 -Sharps +0.0153. 



(3) 



The dispersion of the residuals around the fit is 0.0043, 
which represents only a few percent of the fitted S values. 
Considering intrinsic stellar variability, the quality of the fit is 
probably as good as it can possibly be. As expected, the obtained 
coefficients are very close to a 1:1 relation. Interestingly, the ob- 
served differences between Smw and Sharps (about 0.03) are 
compatible with the estimated scattered light lev el in the H & 
K pas sbands of the Mt Wilson spectrophotometer dDuncan et al.l 
1991). Actually, before correcting the HARPS background, the 
calibration relation was even closer to 1:1, suggesting that both 
instruments roughly have the same amount of scattered light. 
All these results confirm that we have devised an accurate way 
of measuring S mw from HARPS spectra. 

2.4. Conversion to R' HK and metallicity effects 

The Mt Wilson S index provides the Ca II H & K core flux nor- 
malized to the neighboring continuum. This flux contains both a 
photospheric and chromospheric component, whose proportions 
vary as a function of stellar effective temperature and activity 
level. The quantity of interest is usually the chromospheric flux 
alone, which is directly related to the amount of energy that heats 
the chromosphere through the magnetic field. To be able to use- 
fully compare stars between one another, one then has to subtract 
the photospheric component and to normalize the chromospheric 
flux to t he total (bolornetric) luminosity of the star. With this goal 
in mind. lNoves et al.l d!984l) introduced the well-known quantity 
R' UK (and its logarithmic version), which is defined as: 



Ahk = C cf (B - V) ■ S - R phot (B - V), 



(4) 



where C c f (B - V) is a conversion factor that corrects for the 
varying flux in the continuum passbands as a function of B - V 
and normalizes to the bolornetric luminosity, while R v h ot is the 
photospheric contribution in the H & K passbands (also depen- 
dent on B - V). 

The main focus of this paper is to study temporal variations 
in the activity level of our sample stars, and their impact on ra- 
dial velocity measurements and other quantities. In this context, 
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we have to determine which activity-related quantity is the most 
appropriate to consider. We assume that, at the stellar surface, 
the signature of magnetic activity cycles is mainly a variation in 
the average magnetic field strength and filling factor of active 
regions. This is tightly related to the amount of magnetically- 
induced chromospheric heating, measured by R' nK - Further as- 
suming an approximately linear relation between changes in 
chromospheric emission and photospheric line shapes (for small 
variations), we expect that R' UK (and not its logarithm or the S 
index) is the proper quantity to consider. In the remaining of this 
paper we will therefore use R' HK , multiplied by 10 5 for practical 
purposes. The logarithm log/?^ K remains the preferred quantity 
(because in widespread use) when considering only the mean 
activity level of a star. 

We show in Fig.[2]the distribution of mean log R' nK values for 
the 304 stars in the sample, as computed with the above formula. 
A large peak is visible around -4.95. As noted in Sect. 12.11 the 
sample stars have been selected for being slow rotators based 
on v sin i measurements. As a consequence, the sample is biased 
towards old, inactive solar-type stars and cannot be considered 
as representative of the solar neighborhood. As expected, Fig. [2] 
shows a lack of young stars with log R' HK > -4.70, and thus 
does not exhibit the bimodal distribution usually obtained for 
the solar neighborhood (e.g. Vaughan & Preston 1980). As a side 
comment, we note that unbiased planet-search samples such as 
the CORALIE and HARPS volume-limited samples do show a 
bimodal distribution of log R' HK values (Marmier et al., in prep.; 
Lo Curto et al., in prep.). 

Before going further, we examine possible systematics of 
stellar origin in the derivation of R' U ¥ values. The above for- 
mula for R' has been computed by Middelkoop (1982J) and 
iNoves et aD (tl984l) for main-sequence stars only, and taking into 
account a dependence on B - V only. Post-main-sequence evo- 
lution (subgiants) and metallicity are two important factors that 
have an impact on stellar spectra but are not taken into account 
in the R' HK calibration. As already noted by several authors, we 
can therefore expect systematic differences in R' nK values be- 
tween main-sequence stars and subgiants on the one hand, and 
between stars of different metallicities on the other hand. 



It was noted in Sect. 12.11 that 20 stars in the sample are 
actually subgiants based on their position above the main se- 
quence in the H-R diagram. It turns out that all these stars have 
logR' nK < -5.00, with a mean value at -5.11 and a minimum 
value of -5.20. Their distribution is shown as a dashed histogram 
in Fig. |2] We therefore confirm here the well-known, abnor- 
mally low log R' HK values measured on subgiants, caused by their 
higher luminosities and/or lower surface gravities compared to 
main-sequence stars of the same color. We will consider these 
stars separately in the analysis below. 

We now turn to the effects of metallicity. Fig. [3] shows the 
mean R' nK values as a function of stellar metallicity for all stars 
except the 20 subgiants. The lower envelope of the points (de- 
scribing the inactive stars) shows a clear linear trend with metal- 
licity, in the sense that metal-poor stars have higher R' HK val- 
ues than metal-rich stars. This can be understood by considering 
that, at a given B - V color, metal-poor stars are under-luminous 
compared to metal-rich stars. As a consequence, the continuum 
passbands are weaker and the measured S value is larger. This 
leads to a larger R' HK value since the calibration of INoves et alj 
( 1984) assumes that bolometric luminosity (actually o~T^ ff ) only 
depends on B - V. Assuming that most of the metallicity effect 
in R' UK measurements originates in this too simple calibration of 
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values as a function of stellar 



luminosity, we can correct for it by multiplying R' HK with the 
suitable luminosity ratio: 



^H! 



: R u 



T eS (B 



T eS (B 



V,[F£/H])\ 
V,0) J 



(5) 



We now need to express T e ff as a function of B - V and metal- 
licity. Such an expression is actually provided by Sousa et al. 
(2008), obtained from exactly the same data as we are using in 
this paper: 

r eff = 91 14 - 6827 (B-V) + 2638 (B - V) 2 + 368 [Fe/H]. (6) 

The fit is reported to have a standard deviation of the resid- 
uals of only 47 K. Since we are using in this work the same ef- 
fective te mperatures, metalli cities, B - V colors and even stellar 
sample as Sousa et al. (2008), this expression for T e ff is the ideal 
one for our needs. Substituting T & s in Eq.|5] we finally obtain: 



^HK.new -^HK - ' ' l 



368 [Fe/H] 



91 14 - 6827 (B - V) + 2638 (B - VY 



(7) 



The right panel of Fig. |3]is the same figure as the left panel, 
but using the updated calibration. The previously existing trend 
has been successfully corrected, confirming our assumption that 
it was related to an oversimplified B - V vs. r e jf conversion in the 
Noyes et al. calibration. Similarly, we show in the right panel of 
Fig. [2] the corrected distribution of mean logR' nK values for our 
sample. The histogram shows a narrower and higher peak than 
before, centered around -5.00. We interpret this as the sign that 
metallicity differences were artificially smearing the R' HK values 
of stars that actually have very similar activity levels. 

Finally, we plot in Fig. |4]the corrected mean logR' HK val- 
ues as a function of stellar effective temperature. No system- 
atic trend can be seen, confirming again that the R' HK calibra- 
tion takes properly into account the variable stellar properties 
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Fig. 2. Left. Distribution of the mean log R' HK values for the 304 stars in the sample, computed with the Noyes et al. R' HK calibration. 
The shaded area represents subgiants. Right. Same figure, but using the updated R' nK calibration taking into account metallicity. The 
main peak is significantly narrower. 
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Fig. 3. Left. Mean R' HK as a function of stellar metallicity, obtained with the standard Noyes et al. calibration. Metal-poor stars show 
a trend towards higher R' values than metal-rich stars. Right. Corrected mean R' UK - The trend has been successfully suppressed. 



across the mass and metallicity ranges that we consider here, for 
stars on the main sequence. We note however a different dis- 
tribution of logR' HK values between G and K dwarfs: while G 
dwarfs are densely packed around -5.0, K dwarfs (T e ff < 5300 K) 
are clearly more regularly spread over the range -4.7 to -5.0. 
This can possibly be interpreted as the signature of a slower de- 
crease in activity level with age among K dwarfs, as discussed by 
Mamajek & Hillenbrand (2008). These authors provide, in their 



Fig. 11, a calibration of the log^ K evolution as a function of 
age and spectral type. Assuming that G and K dwarfs in the so- 
lar neighborhood have roughly the same age distribution (~2-8 
Gyr for our sample), we expect that most G dwarfs will have 
already reached the logR' HK = -5.0 level, while K dwarfs are 
still spread between -4.7 and -5.0. The observed accumulation 
of G dwarfs at -5.0 suggests however that stars do not smoothly 
evolve towards even lower activity levels once they have reached 
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this floor, but rather remain at -5.0 for an extended period of 
time. 



2.5. Radial velocities and line shape parameters 

The high-precision radial velocities and line shape parameters 
used in this paper are those obtained with the standard HARPS 
pipeline (see Lovis et al., in prep., for a more detailed descrip- 
tion). They are all derived from the cross-correlation function 
(CCF) that is built using a binary template containing several 
thousands of spectral lines. Radial velocities are then obtained 
by simply fitting a Gaussian to the CCF profile and measuring 
its centroid. The fit also yields the FWHM of the CCF and its 
contrast, which we use here as spectroscopic indicators of stel- 
lar activity. Moreov er, we also compute the CCF bisector span 
( Ouel oz et al.l2 00 1 ), which represents the velocity difference be- 
tween the upper and lower parts of the CCF, and is useful to 
characterize line shape asymmetries. 

HARPS not only delivers precise radial velocities at the 
level of ~1 ms _I and better, but also spectroscopic indicators 
of the same quality. This is extremely valuable in the context 
of exoplanet searches since it allows us to monitor stellar activ- 
ity simultaneously with the radial velocity measurements, and 
this not only for active stars but also for 'quiet' ones with 
log R' nK around -5.0. Several recent papers have already high- 
lighted the us e that can be made of th ese high -precision indi- 
cators (Oueloz et a l]l2q09HSantos et alJl2010t Lovis et ai]l201lt 
iDumusque et alj|201 latlPepdl201 lh . 



3. Statistics of Ca II H & K variability 

3. 1 . Raw dispersions and measurement errors 

Before searching for signals of stellar origin in chromospheric 
emission measurements, we first examine the raw dispersions of 
the Ri, K time series and try to characterize instrumental errors. 
Fig. |5]shows the distribution of the R' HK standard deviations for 
the 304 stars in the sample. The peak of the distribution is lo- 
cated at ~0.02. A tail of higher- variability stars is present, either 
due to short-term variability or magnetic cycles as we will see 
below. On the other side of the distribution, we note that several 
well-measured stars exhibit a scatter around 0.01. As an exam- 
ple, Fig.|6]shows the time series of the well-known standard star 
t Ceti (HD 10700), which comprises 153 data points spread over 
more than 7 years. The dispersion of the data is only 0.0089 in 
R' nK , corresponding to a relative variation around the mean of 
0.93%, and to a relative variation of the S index of 0.35%. This 
demonstrates both the high instrumental precision reached by 
HARPS and the remarkable stability of this star. Clearly, t Ceti 
does not show any magnetic activity cycle over the time span 
of the observations, confirming its previously reported status of 
very quiet star. 

The determination of reliable error bars on the R' HK measure- 
ments is a difficult issue. On the one hand, one has to disentangle 
instrumental errors from intrinsic stellar variability and on the 
other hand, instrumental errors may be dominated by systematic 
effects that are difficult to estimate. Because the SNR in the Ca II 
H & K line cores is usually low (less than 5:1 per extracted pixel 
for the fainter stars in the sample), and because of the large SNR 
contrast with the continuum, S-index measurements are sensitive 
to several instrumental effects like non-perfect diffuse light sub- 
traction, residual cosmic hits, CCD defects, CCD charge transfer 
inefficiency and non-linearity, etc. (see Sect. l2.21 >. 
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Fig. 5. Distribution of the R' nK dispersions for the 304 stars of 
the sample. 
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Fig. 6. Time series of R' HK measurements for the standard star r 
Ceti (HD 10700). 



The observed R' HK dispersions are typically 2-4 times larger 
than the estimated photon noise levels for the sample stars, most 
of which are between 0.002 and 0.020 (median value 0.007). 
Plotting the normalized standard deviations as a function of pho- 
ton errors, we note that the lower envelope of the points departs 
from 1.0 when photon errors become smaller than ~0.007. This 
suggests that either all stars are intrinsically variable at least at 
that level, or that we have hit an instrumental noise floor. Given 
the various effects mentioned above, we are more inclined to be- 
lieve in the instrumental explanation. Consequently, we quadrat - 
ically add an additional noise term of 0.007 to all error bars on 
R' nK in the following analysis (we note that the precise value of 
this noise floor is unimportant for the magnetic cycle search). 
Towards lower SNR, the photon errors seem to correctly de- 
scribe the lower envelope of standard deviations, at least up to 
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a photon noise level of -0.03 (the worst case for our sample). It 
is likely that other effects come into play at even lower SNR (e.g. 
non-perfect background subtraction), but that regime is beyond 
the scope of this paper. 

In conclusion, we reach a long-term precision of ~0.01 in 
R' HK on bright stars, corresponding to a relative precision os /S 
of ~0.35%. Interestingly, these HARPS measurements are about 
a factor of 3 more precise than those of the Mt Wilson H & 
K project, which rea ch a long-term precision of about 1 % in S 
(Bali unas et al.lfl995l) . To the best of our knowledge, the long- 
term chromospheric data presented here are likely the most pre- 
cise ever reported for stars other than the Sun. As a final com- 
ment on measurement precision, we note that in the case of the 
Sun, typical, small active regions induce a modulation in chro- 
mospheric emission of ~l% as they rotate with the star. The pre- 
cision reached by HARPS should therefore allow us to be sen- 
sitive to individual active regions crossing the disk of solar-type 
stars, provided they are observed at a suitable cadence. Many ro- 
tational modulations are indeed readily visible in HARPS R' nK 
data, which will be the topic of a subsequent paper. 



3.2. Method and classification 

The large dataset at our disposal allows us to perform a global 
study of the behavior of chromospheric emission in old solar- 
type stars as a function of time. In this exercise, the main diffi- 
culty arises from the extremely inhomogeneous time series for 
individual stars, i.e. the varying number of observations per star 
and the time span of the observations. Indeed, observations were 
scheduled as a function of the radial velocity behavior of each 
star, with planet-host stars receiving significantly more attention 
than RV-stable stars. A global analysis of chromospheric vari- 
ability has therefore to be robust enough to deal with scarce as 
well as dense time series in a consistent way. 

As a first step, our main goal is to classify stars into broad 
categories according to their Ca II H & K variability. We char- 
acterize variability in two different ways: 1) a search for signif- 
icant long-term periodicities in the time series, and 2) a short- 
term variability criterion based on the standard deviation of the 
measurements. 
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Fig. 7. Top: GLS periodogram of the R' HK data for the star 
HD 7199. The dashed line denotes the 1% false-alarm probabil- 
ity threshold. Bottom: Time series of the R' HK data with the fitted 
magnetic cycle model overplotted. 



3.2.1. Periodicities 

To search for period icities, we use the generalized Lomb-Scargle 
periodogram (GLS, IZechmeister & Kursterll2009l) and consider 
a given periodogram peak as significant when its power exceeds 
the 1 % false-alarm probability (FAP) threshold. FAP levels are 
computed by performing random permutations of the data while 
keeping times of observations fixed. This procedure provides 
a robust way of finding significant periodicities in a very in- 
homogeneous dataset in terms of number of observations per 
star. We fit magnetic cycles with a simple sinusoid at the period 
of maximum power, as provided by the GLS. As an example, 
Fig. |7] shows the periodogram and fitted time series for the star 
HD 7199, which exhibits a strong magnetic cycle at a period of 
2760 days (see Dumusque et al. 2011a, for a detailed study of 
this star). 

We further split stars with a detected signal into two groups: 
short-term and long-term periodicities, with a limit between the 
two at P = 2 yr. We choose 2 yr for two main reasons: 1) it is an 
intermediate timescale separating the stellar rotation timescale 
on the one hand and stellar magnetic cycles on the other hand, 
and 2) it minimizes the risk of misclassifying stars because of 



aliasing problems related to various aliases in the window func- 
tion of the observations. Putting the limit at shorter periods cre- 
ates a risk of misclassifying stars that show strong rotationally- 
induced variability and poor sampling, while putting the limit at 
longer periods may miss true activity cycles (i.e. not related to 
rotational modulations) that exist at timescales of a few years. 
In brief, poor sampling often prevents us from reliably detect- 
ing magnetic cycles below ~2 yr (if they at all exist) because 
it is then difficult to disentangle them from rotationally-induced 
short-term variability. We therefore focus hereafter on magnetic 
cycles in the 2-11 yr period range. 

For all stars that exhibit no significant long-term modula- 
tions, we compute detection limits on a possible sinusoidal sig- 
nal at long periods. This allows us to use all available infor- 
mation in our dataset, i.e. to put constraints on the existence of 
magnetic cycles for apparently constant stars and stars with only 
few measurements. The output of these computations is an upper 
limit on the R' nK semi-amplitude that cannot be excluded by the 
data. In order to obtain a unique limit for all long-period signals, 
we average the upper limits over periods of 2-1 1 yr. We choose 
1 1 yr because the typical time span of the observations (~6-7 
yr) allows us to have meaningful upper limits roughly up to this 
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Fig. 8. Detection limits for sinusoidal signals in R' nK measure- 
ments for the standard star t Ceti (HD 10700). The step-like 
curves are due to the step size of 0.001 used in the computations. 



timescale, and because it is the period of the solar magnetic cycle 
(to be able to compare with the Sun). The Sun exhibits a roughly 
sinusoidal variation in R' HK with a semi-amplitude A of -0.26. 
As we will see below, the sample stars show magnetic cycles 
with a variety of semi-amplitudes, including very low ones. It 
is therefore difficult to define a threshold amplitude in R' nK that 
would exclude the existence of a magnetic cycle. In this work we 
arbitrarily choose A - 0.04 as the minimum semi-amplitude for 
a Targe-amplitude' magnetic cycle. This corresponds to about 
15% of the solar R' HK semi-amplitude. This choice is motivated 
by the average sensitivity of our survey: a much lower thresh- 
old would reject too many stars and the global results would be 
strongly affected by small-number statistics. 

In computing detection limits, an artificial signal is injected 
into the data and its detectability is checked. There are three 
free parameters in this exercise: the period, amplitude and phase 
of the signal. While the goal is to derive a limiting amplitude 
as a function of period, the phase is an extra parameter that 
must be probed independently, adding an extra dimension to the 
problem: detection limits should actually be computed for every 
phase value. At each period, one can then derive the distribution 
of detection limits for all phases and, assuming all phases are 
equiprobable, finally assign a probability of detecting a given 
amplitude at a given period. While this treatment of phase may 
not be necessary for well-sampled time series with no signifi- 
cant gaps, it becomes important for our sample of irregularly- 
observed stars because detection limits at a given period often 
vary widely as a function of phase. In summary, for stars with no 
significant long-term signals, we compute detection limits for a 
set of 12 equidistant phases at all periods. From this, and given 
a semi-amplitude Aq, we can compute a lower limit to the prob- 
ability that no cycle of semi-amplitude Aq or larger exists: 



Prob (A < Aq) > 



2n Ja,„ 



d<p, 



(8) 



where A\i m (<f>) is the detection limit as a function of phase </>, 
which takes values between and 27r. 

Fig. |8] shows the results of the detection limits calculations 
for r Ceti. The different curves denote the 12 equidistant phases 
simulated at each period. The upper curve represents the detec- 
tion limit for the 'worst' phase, i.e. the minimum semi-amplitude 



that can be detected whatever the phase of the signal. The lower 
curve represents the minimal semi-amplitude that can be de- 
tected for the most favorable case, i.e. for the best 8% (=1/12) of 
all phases. At a period of 15 yr, we see that a semi-amplitude of 
0.01 can be excluded for many phases, but not for all. Counting 
how many curves are crossed from semi-amplitude to 0.01, 
we can estimate that a signal with semi-amplitude 0.01 at 15 yr 
can be excluded with a probability larger than 67% (=8/12). As 
another example, we immediately see from Fig. [8] that a semi- 
amplitude of 0.04 (the limit for a Targe-amplitude' magnetic cy- 
cle) can be fully excluded for periods between 2-20 yr. 

In summary, we define the following variability classes re- 
garding magnetic cycles: 

- Stars showing a magnetic cycle, i.e. significant long-term 
(P > 2 yr) modulation; among these, Targe-amplitude' mag- 
netic cycles are those with R' HK semi-amplitude A > 0.04. 

- Stars without a magnetic cycle: the detection limit on the 
R' nK semi-amplitude is < 0.04 for periods between 2-11 yr. 



3.2.2. Short-term variability 

In this work we mainly focus on long-term modulations in chro- 
mospheric emission, i.e. magnetic cycles. Therefore, stars ex- 
hibiting significant long-term periodicities in R' HK will receive 
particular attention in the following sections. However, signif- 
icant variability also occurs on the short term, and we attempt 
to quantify it through the standard deviation of the measure- 
ments, computed after subtracting the best-fitting long-term sig- 
nal (P > 2 yr), whether significant or not. Short-term variability 
in R' nK contains a wealth of information about stellar rotation 
that deserves a separate study. Again, the inhomogeneity of the 
time series makes a global analysis more complicated, since the 
cadence of the observations plays a major role in the ability to 
characterize rotational modulations. Whereas the best-observed 
stars may reveal their rotation periods by a simple look at the 
R' HK curve, stars with few measurements will only show scatter 
with undetermined characteristic timescale. In this context, the 
(short-term) standard deviation of the measurements provides a 
useful diagnostic that is not very sensitive to the number of data 
points per star. 



3.3. Results 

The global results of the R' HK analysis are given in Tables [2] and 
|3]for all stars in the sample. Table |2]contains the stars showing a 
significant magnetic cycle. We fit magnetic cycles with a simple 
sinusoid at the peri od of maximum power, as pr ovided by the 
GLS periodogram (Zechmeister & Kurster 2009). Confidence 
intervals for the period are obtained from Monte Carlo simu- 
lations: in each trial, different noise realizations on individual 
data points are drawn from a Gaussian distribution with a stan- 
dard deviation equal to the measured short-term scatter in the 
time series, and the GLS periodogram is then recomputed on the 
modified dataset. The distribution of periods of maximum power 
yields the desired confidence intervals for the period. 

Provided for each star in Table [2] are: the number of obser- 
vations, their time span, the mean and median log R' nK values, 
the raw and short-term R' UK standard deviations, the magnetic 
cycle period with 68% confidence intervals, the mean log R' HK 
value y as obtained from the magnetic cycle fit, and the cycle 
semi-amplitude A. Finally, estimates of the stellar rotation pe- 
riod and age are also provided, based on the activity-rotation-age 
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calibration of Mamajek & Hillenbrand (2008) and using y as the 
long-term mean logR' HK value. 

Table [3] provides results for stars with no detected magnetic 
cycle. The first 7 columns of the table are the same as for cycling 
stars. Then the detection limit on the R' nK semi-amplitude An m 
is provided, averaged over periods of 2-11 yr (maximum value 
over all phases). The following column gives the minimum prob- 
ability that no cycle exists, taking into account the variable sen- 
sitivity as a function of phase (see Eq.[8]). Finally, estimated rota- 
tion periods and ages are also given, based on the mean log R' HK 
value. 



3.4. Magnetic cycles 

3.4.1. Occurrence 

We now examine the results of the search for magnetic cycles. 
Among 284 main-sequence FGK stars, we find 99 stars with 
a detected magnetic cycle, and 20 stars where a magnetic cy- 
cle with A > 0.04 over 2-11 yr can be fully excluded. This 
means there are 165 stars with no detected cycle, but for which 
the existence of a cycle cannot be completely excluded. Among 
the 99 stars with a magnetic cycle, 16 of them have R' HK semi- 
amplitude smaller than 0.04. Therefore, there are 83 stars with a 
large-amplitude cycle and 36 stars with no large-amplitude cy- 
cle. From this we deduce that about 70% (83/119) of the sam- 
ple stars exhibit a large-amplitude cycle, while 30% (36/1 19) do 
not. These estimates assume that there is no observational bias 
related to the classification of stars in either category, e.g. that 
stars in either category have not been preferentially observed. 
Observations have been primarily scheduled as a function of the 
radial velocity variability of stars, and only indirectly according 
to their activity behavior. We however need to be cautious about 
observational biases because it is generally easier to detect an ex- 
isting cycle than to completely exclude any cycle with amplitude 
A > 0.04. Indeed, the latter case requires more data points on av- 
erage than the former one because of the need to put constraints 
at all phases, and the increased cycle detectability as amplitude 
increases. The fraction of 30% of non-cycling stars is therefore 
probably a lower bound. More conservatively, one should first 
provide the occurrence rates normalized over the whole sample: 
29% (83/284) of the sample stars have a large-amplitude cycle, 
13% (36/284) do not, and 58% (165/284) have an undetermined 
status. 

We try to further refine these statistics of occurrence by using 
the information from as many stars as possible, i.e. taking into 
account the detection limits on R' UK signals also for stars where a 
magnetic cycle larger than 0.04 cannot be fully excluded (i.e. at 
all phases). We follow the methodology described in Sect. 13.2. Tl 
and compute for each sample star the lower limit to the proba- 
bility that no cycle larger than 0.04 exists for periods between 
2-11 yr (Eq.[8]l. By summing all these probabilities one can de- 
rive another estimate of the proportion of quiet stars in the whole 
sample. We obtain that at least 26% of sample stars do not have 
a magnetic cycle, to be compared to at least 29% which do have 
one (and 45% remain unclassified). In this case, the 'classifica- 
tion bias' has shifted in the other direction since we used stars 
with less data points for the non-cycling category, and the 26-to- 
29% ratio can be considered as an upper bound to the fraction 
of non-cycling stars. We conclude that between 30 and 47% of 
older solar-type stars in the solar neighborhood do not have any 
magnetic cycle with A > 0.04, while between 53 and 70% do 
have one. 
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Fig. 9. Distribution of magnetic cycle periods for the 99 stars ex- 
hibiting significant long-term modulations in R' UK - Periods are 
poorly constrained beyond ~4000 days due to the limited du- 
ration of the survey. Unconstrained long-term trends are repre- 
sented in the last bin at 8000 days. 
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Fig. 10. Distribution of magnetic cycle semi-amplitudes for the 
99 stars exhibiting significant long-term modulations in R' HK - 



3.4.2. Periods and amplitudes 

Fig- and Fig. [TOl show the distributions of magnetic cycle pe- 
riods and semi-amplitudes, respectively. Perhaps surprisingly, 
the period distribution shows a peak at relatively short periods 
(~3-5 yr) and a sharp decrease at long periods. Our survey is 
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not long enough to properly constrain periods beyond ~4000 
days, but long-term trends are nonetheless easily detected if they 
are of sufficient amplitude. These are shown in the last bin in 
Fig. [9] Even if all these unconstrained long-term trends were 
better characterized, they would not be in sufficient numbers to 
significantly change the overall shape of the period distribution. 
The Sun, with its period of 11 yr (4000 d) appears to have a 
rather long cycle compared to other solar-type stars (80% of the 
stars have a period shorter than the Sun). As for cycle semi- 
amplitudes (Fig. [TOt. we note a broad peak at ~0.1 in R' HK and 
a large number of low-amplitude cycles. Again, the Sun appears 
to have a rather large A of 0.26 compared to other solar-type 
stars (82% have a lower semi-amplitude). Finally, we note that 
the first bin in Fig. [10] contains semi-amplitudes smaller than 
0.04, i.e. cycles that are not considered Targe-amplitude' in the 
present work. We nevertheless see that several stars show vari- 
ations at that level, but these cycles are difficult to measure and 
require a large number of data points to be detected. 



3.4.3. Correlations with activity level and stellar properties 

We further searched for possible correlations between cycle pe- 
riod, amplitude, mean activity level and effective temperature. 
We found nothing obvious between period and amplitude, ex- 
cept perhaps a lack of large-amplitude cycles at short periods. 
When examining period vs. mean activity level, it seems that 
stars more active than log R' nK - -4.75 only exhibit short periods 
(P < 3000 d), although this trend remains uncertain because of 
small-number statistics. Overall, it is difficult to analyze cycle 
periods in more details because of the rather large error bars af- 
fecting many periods. A clearer picture emerges when plotting 
cycle amplitude as a function of mean activity level (see Fig.fTTI). 
Inactive stars with logR' nK = -5.0 only show low cycle ampli- 
tudes, while more active stars tend to show larger amplitudes on 
average, although there is a rather large scatter around this trend. 
We also plot in Fig. [T2l cycle semi-amplitude vs. effective tem- 
perature. We see that the largest amplitudes are concentrated in 
the range 5000-5600 K, with a decreasing trend on both sides 
of the temperature scale. In particular, hotter stars do not show 
semi-amplitudes larger than ~0.2. 

We now examine how stars with and without a magnetic cy- 
cle are distributed in a H-R diagram and in a T e ff vs. logR' HK 
plot. Fig. [J_3] shows a H-R diagram of the 284 main-sequence 
stars in the sample, with a different symbol used for a) stars with 
a large-amplitude cycle, b) stars with no large-amplitude cycle, 
and c) stars with insufficient observations to be classified. We 
immediately see that almost all non-cycling stars are found at 
higher temperatures (T^ > 5300 K), while stars with large am- 
plitudes seem to exist at all temperatures. In the T e ff vs. log R' nK 
plot (Fig. [141. we see that all non-cycling stars are concentrated 
in the low-activity, high-temperature regime. In particular, there 
is no constant star above logR' nK = -4.94. Stars with a large- 
amplitude cycle are distributed roughly homogeneously in tem- 
perature at larger mean activity levels, although there are some 
cycling stars below logR' HK = -5.0. 

We can think of two possible interpretations for these results. 
Firstly, considering the mean log R' HK value as an age indicator 
for solar-type stars, these findings would indicate an evolution 
in magnetic cycles characteristics as stars age, in the sense of a 
decrease in cycle amplitude with age. This interpretation is also 
supported by Fig.Qj] and the existence of many low-amplitude 
cycles (Fig. ITOl , This picture fits well with the known paucity 
of magnetic cycles in subgiants (see Sect. 13.4.41 ). However, a 
second explanation might also be possible: since we know that 
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Fig. 11. Magnetic cycle semi-amplitude as a function of mean 
activity level. 
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Fig. 12. Magnetic cycle semi-amplitude as a function of stellar 
effective temperature. 



the Sun undergoes periods of extended activity minima (e.g. the 
Maunder minimum), it is possible that the non-cycling stars in 
our sample are in fact in a temporary state of minimum activity. 
This is supported by the fact that some stars have been caught in 
a tran sition state between cycling and non-cycling behaviors (see 
e.g. Baliu nas et al.lll995l) . In the end, both explanations may be 
simultaneously true, and stars may actually go through a phase 
of transient magnetic cycles as they age, before these cycles 
completely disappear when leaving the main sequence. 
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Fig. 13. H-R diagram showing stars with a large-amplitude mag- 
netic cycle (red), stars with no large-amplitude cycle (blue), and 
stars of undetermined status (black). 



Fig. 14. T e g vs. logR' HK plot showing stars with a large- 
amplitude magnetic cycle (red), stars with no large-amplitude 
cycle (blue), and stars of undetermined status (black). The posi- 
tion of the Sun is shown as a G symbol. 



It must be noted that this discussion mainly concerns G 
dwarfs, since we find only one non-cycling star among K dwarfs 
(T e ff < 5300 K). We have already noted that the distribution of 
mean log R' MK values varies markedly between G and K dwarfs 
(see Sect. 12.41 ). Whereas log R' UK values show a clear peak at -5.0 
for G dwarfs, they are much more spread between -4.7 and -5.0 
for K dwarfs. The lack of non-cycling K dwarfs may be due to 
the fact that these stars are on average not old enough to have 
reached the logR' HK a; -5.0 level, which may be required to halt 
magnetic cycles. 



3.4.4. Subgiants 

It has long been known that stars evolving off the main sequence 
tend to show part icularly low levels of chromospheric activity 
(e.g. Wright 2004), and also less long-term variability than main- 
sequence stars. We confirm these findings here. As noted earlier, 
all 20 subgiants in the sample have a mean logR' nK value be- 
low -5.0. More importantly, only one of them shows a magnetic 
cycle (semi-amplitude A = 0.055), while 5 stars have no cycle 
at all (A < 0.04). Moreover, many of the 14 stars not classified 
in either category have high probabilities to have no magnetic 
cycle. Summing these probabilities yields a global proportion 
of between 83-93% of non-cycling stars among subgiants, to be 
compared to 30-47% for main-sequence stars. We therefore con- 
clude that most subgiants have indeed reached a stage of long- 
term stability in their chromospheric emission. 



3.4.5. The case of the Sun 

We can now try to consider the Sun within the broader context of 
our sample of solar-type stars. We show in Fig.[14]the position 
of the Sun at r eff = 5770 K and logR' HK = -4.91. It is mainly 
surrounded by cycling stars in this figure, but is actually not far 
from the boundary between mainly cycling and mainly constant 



G dwarfs, located at logR' HK w -4.95. One could therefore con- 
jecture that the Sun will soon transition to a permanently quiet 
state, perhaps within ~l-2 Gyr. On the other hand, the present- 
day solar cycle semi-amplitude of 0.26 appears rather large com- 
pared to other stars at the same activity level and effective tem- 
perature (see Fig.[TT1and[T2l. However, one should keep in mind 
that there are significant variations in amplitude from one solar 
cycle to the next (up to 50% and more), and that solar activity 
has been rather high since 1950 compared to previous periods, as 
indicated by sunspot data. On average, the Sun is therefore likely 
to have a cycle amplitude that is typical of stars in its activity and 
temperature range. It is also conceivable that the Sun is presently 
in an unstable state with activity cycles that are highly variable 
in amplitude, and that this could be the prelude to a transition to 
a permanently quiet state. 



3.5. Short-term variability 

Besides the long-term variations in chromospheric emission 
caused by magnetic cycles, solar-type stars also show (some- 
times prominent) short-term variability caused by active regions 
rotating with the star. For all stars in the sample, we characterize 
this short-term variability by computing the R' HK standard de- 
viation, after subtracting the best-fitting long-period signal with 
P > 2 yr (whether significant or not). We then search for cor- 
relations between R' HK variability and stellar parameters. The 
only clear dependence we could find is with the mean activ- 
ity level of the stars. Fig. Q3] shows the dispersion in R' HK as 
a function of logR' HK itself. As can be seen, there is a dense 
clump of low-variability stars with log R' HK between -5.1 and - 
4.9, and a marked trend towards high variability at higher ac- 
tivity levels. However, there are many variable stars at logR' nK 
- -4.90 already. Actually, the upper envelope of the points in- 
creases sharply from -5.00 towards higher activity levels. The 
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Fig. 15. R' HK short-term dispersion as a function of mean activity 
level log R' HK . 

lower envelope appears flat until about -4.95, where it starts in- 
creasing steadily with increasing activity level. The flattening 
at low activity may be due either to a genuine stellar variabil- 
ity floor or, more probably, to residual instrumental systematics. 
Incidentally, we note that the Sun reaches an activity level of 
about -4.97 at cycle minimum, when almost no sunspots are vis- 
ible at the solar surface. It may be that this level represents a 
universal threshold that marks the appearance of significant ac- 
tive regions on solar- type stars. 

We also note that the scatter around the overall trend in 
Fig.[T"5lmav be caused at least partly by inclination effects, i.e. by 
the different space orientations of stellar rotation axes. Indeed, 
stars seen pole-on are expected to show less rotationally-induced 
variability since active regions maintain the same projected area 
as seen from the observer's standpoint. 

In summary, it appears that the mean activity level, expressed 
by logR' nK , is a good indicator of magnetically-induced variabil- 
ity on the short-term as well as for long-term activity cycles, 
since these also tend to increase in amplitude with the mean ac- 
tivity level. 

In the context of radial velocity searches for exoplanets, 
more in-depth studies of short-term variability are fundamental 
to understand radial velocity jitter induced by stellar activity. It 
is clear that the degree of R' HK short-term variability is linked to 
the degree of RV, FWHM, contrast and BIS variability. However, 
the exact relationships between all these quantities on the short 
term are complex and beyond the scope of this paper. Several 
well-sampled time series are available among the stars in this 
sample and deserve individual investigations. These will be the 
topic of a subsequent paper (Dumusque et al., in prep.). 

4. Impact of magnetic cycles on radial velocities 
and line shape parameters 

We now want to study the possible impact of magnetic cycles on 
radial velocities, CCF FWHM, CCF contrast, and CCF bisec- 
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Fig. 16. Time series of R' nK , RV, FWHM, contrast and BIS mea- 
surements for the star HD 21693, showing clear correlations be- 
tween all quantities. 



tor span. We consider the 99 main-sequence stars that have been 
found to have a magnetic cycle and try to find correlations be- 
tween R' UK and the other quantities. In several cases, correlations 
are actually obvious to the eye. An example is shown in Fig. [TBI 
where the time series of the five quantities are displayed for the 
G8V star HD 21693. As can be seen, the RV, FWHM, contrast 
and BIS measurements clearly vary in phase with R' nK - Similar 
behaviors are apparent on many other stars, but often less clearly 
so because of poorer sampling. We therefore need robust meth- 
ods to study correlations in the whole sample of stars. Two of 
them are described below. 



4.1. Blind fitting of the R' HK -dehved model 

One possible method is to fix the period and the phase of the pu- 
tative RV, FWHM, contrast or BIS signal to the values found 
when fitting the magnetic cycle, and then compute via linear 
least-squares the best-fitting mean levels and semi-amplitudes 
on the RV, FWHM, contrast and BIS data. As for magnetic cy- 
cles, this assumes that we are looking for sinusoidal variations 
in these data, and that these variations are in phase with the R' nK 
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data. The advantage of this method is its robustness, since only 
two free parameters are fitted on the whole time series, and the 
long-period sinusoidal model is relatively insensitive to short- 
term scatter related to rotational modulations. In this way, effects 
from magnetic cycles on the one hand, and rotational modula- 
tions on the other hand can often be distinguished even if sam- 
pling is rather poor. 

Following this scheme, for each star with a magnetic cycle 
we obtain a best-fit mean value and semi-amplitude for the RV, 
FWHM, contrast and BIS data, at the period and phase given by 
the R' UK data. Then the sensitivity of each of these quantities to 
magnetic cycles is obtained by simply dividing the RV, FWHM, 
contrast and BIS semi-amplitudes by the R' nK semi-amplitude. 
Uncertainties in the sensitivities are derived from the covari- 
ance matrix of the least-squares fits. In this process we artifi- 
cially inflate error bars on individual measurements to bring the 
reduced^- 2 to 1, so as to "model" the short-term scatter in each 
quantity as a random Gaussian noise source. This approximation 
may be risky on a case-by-case basis since rotational modula- 
tions induce signals that are time-correlated and non-Gaussian. 
However, these aspects have only little effect in our case since 
we are fitting only the mean level and amplitude of a long-period 
signal, and we are mainly interested in the global results over the 
whole sample rather than to individual stars. 

Fig. [TTl shows the obtained sensitivities of the RV, FWHM, 
contrast and BIS, plotted as a function of effective temperature. 
We emphasize that all stars are included in the procedure, with 
no attempt to filter out objects that show variability manifestly 
unrelated to magnetic cycles, including RV signals caused by or- 
biting companions, short-term activity or instrumental systemat- 
ics. In spite of this, a clear dependence on T e g can be seen for 
all quantities. Besides T e ff, the other fundamental stellar param- 
eter that is likely to play a role is metallicity. We proceed to fit 
linear models in T e ff and [Fe/H] to the sensitivities. Models com- 
bining a dependence on both parameters are found to yield the 
best reduced^- 2 , and second-order polynomials are found to be 
sufficient for both parameters (going beyond would increase the 
reduced^ 2 ). We finally obtain the following empirical models: 



C, 



RV 



(19.04 ±0.62) + (3.20 + 0.20) 

+ (4.5 + 3.1) 

+ (24.3 + 5.2) • [Fe/H] 2 



10~ 6 • f 2 ff + (17.8 



C 



FVV 



(83.2 + 1.4) - (6.46 + 0.36) • 10~ 2 • T eff 

- (3.71 + 0.79) • 10~ 5 • 1 2 ff + (33.7 ± 6.8) • [Fe/H] 

+ (38 ± 17) • [Fe/H] 2 (10) 



Ceo = (-0.703 + 0.012) + (3.9 + 3.4) -10 



• T e fi 



+ (5.28 : 

-(0.15: 



1-7 *2 



0.67) • 10"' • T e z ff + (0.333 ± 0.055) • [Fe/H] 



0.13) -[Fe/H] 2 



Cbis = (20.18 + 0.44) + (7.3 ± 1.3) ■ 10~ 3 • f eff 

- (7.8 + 2.4) ■ 10~ 6 ■ T 2 ff + (8.2 + 2.3) • [Fe/H] 
+ (13.7 + 5.1) -[Fe/H] 2 



(11) 



(12) 



In these equations, Crv, Cfw, Cco and Cbis are the sensitiv- 
ities expressed in m s (or percent for the contrast) per unit of 
R' nK , and f e jf = r e ff-5400K (chosen to minimize covariances 
between polynomial coefficients). Fig. [17] shows these models 
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Fig. 17. Sensitivity of RV, FWHM, contrast and BIS measure- 
ments to magnetic cycles (expressed by R' nK variations), as a 
function of T e g. 
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overplotted on the data for zero metallicity. Some of the scatter 
seen in the figure is in fact due to the effects of metallicity. 

A dependence on stellar properties, in particular T e ff and 
[Fe/H], is expected since photospheric lines will respond differ- 
ently to magnetic cycles depending on physical conditions and 
the details of the atmospheric structure. An important issue aris- 
ing here is whether the simple dependencies derived above are 
able to capture the full effects of magnetic cycles on spectral 
lines, or whether other factors play an important role as well. A 
way to estimate this is to look at the reduced x 2 of the best-fit 
models to check if some variability remains unaccounted for. A 
related question is whether we should filter out some stars before 
performing the fit if they show additional variability unrelated to 
magnetic cycles. Taking all stars into account, we obtain reduced 
X 2 values of 6.79, 2.51, 5.64 and 1.91 for RV, FWHM, contrast 
and BIS, respectively. Clearly, additional variability is present 
in radial velocities and contrast, while FWHM and BIS data are 
closer to a reasonably good fit. The RV result is not surprising 
since we are considering here raw measurements, which in par- 
ticular still contain the signals of many planets. We have no ex- 
planation for the rather poor^ 2 value for the contrast, except that 
this parameter is the most sensitive to instrumental systematics. 

We examine below the possibility to remove the RV signals 
of known planets before deriving RV sensitivities. But as a first 
step, we try to estimate the impact these signals might have on 
our fitted sensitivity curve (Eq.|9]). To do this, we artificially re- 
move from the data the three most significant outliers as well 
as all stars with R' HK semi-amplitudes smaller than 0.04 (which 
generally have less precise sensitivities), and we refit this filtered 
dataset. The obtained coefficients are indistinguishable from the 
original ones in Eq. [9] This shows that the exact choice of stars 
to include in the fit, and the presence of planet-induced RV sig- 
nals, have a negligible effect on the final RV sensitivity curve. It 
is mainly due to the fact that error bars on individual sensitivi- 
ties usually take into account the mismatch between RV and R' HK 
data if signals unrelated to magnetic cycles are present. In other 
words, a significant outlier can only arise if a planet induces a 
strong RV signal exactly in phase with the magnetic cycle of its 
host star, which is obviously quite unlikely. 

We conclude that the obtained sensitivity curves in Eq.|9]are 
valid, even if they were derived based on uncorrected RV data. 
The relatively good reduced x 2 values for the FWHM and BIS 
data also indicate that the simple r e ff and [Fe/H] dependence is 
able to capture most of the variability seen in the sensitivities, 
and that other stellar properties probably have a much smaller 
effect. 

4.2. Data binning and subtraction of planetary signals 

In parallel, we also try an alternative approach for the RV sensi- 
tivities. We search for correlations directly through linear regres- 
sion between RV and R' nK data, using the Pearson correlation 
coefficient to estimate the quality of the correlation. Before that 
however, considering that RV data are affected by the presence 
of orbiting companions, we first remove the signals of planetary 
origin from the RV data to obtain a more reliable RV - R' nK cor- 
relation. We also restrict ourselves to stars showing a magnetic 
cycle with A > 0.04 (83 stars out of 99). Moreover, to have suffi- 
cient sampling and minimize the impact of short-term variability 
we further refine the selection as follows: 

- Stars must have more than 3 years of measurements 

- Measurements are binned over 3 months, with a required 
minimum of at least 3 measurements per bin 



- A total of at least 5 bins must be available 

After this selection we are left with 44 stars, which all show 
clear magnetic cycles and offer a rather good sampling over more 
than 3 years. The following step consists of subtracting the RV 
signals of known planets, without erasing a possible signal com- 
ing from magnetic cycles. For the already known planets we use 
the orbital parameters found in the literature. There are also sev- 
eral new objects in this HAR PS sample, for which the publica- 
tion is in preparation (|Pepell201lt IDumusau e et al.ll201 lai Udry 
et al., in prep; Queloz et al., in prep.). In these cases we take 
the preliminary orbital elements from these unpublished works. 
Usually, the magnetic cycle is first fitted with a sinusoid on the 
activity index, and the obtained model is then fitted to the RV 
data together with the Keplerians, just leaving the amplitude as 
a free parameter (see iDumusque et al.ll201 lal for an example). 
The RV - R' UK correlation slope is then obtained using the RV 
residuals (i.e. all planets subtracted). We plot in Fig.[18]the RV 

- R' nK correlation slope as a function of effective temperature for 
the 44 selected stars. 

As can be seen in Fig. [18] the RV - R' nK correlation slope is 
well correlated to effective temperature. However, some stars fall 
quite far from the fitted quadratic model. Most stars which seem 
to be outliers happen to have a square symbol in Fig. [18] which 
represents stars with a maximum activity value higher than 
logR' HK - -4.75. This corresponds to activ e stars according to 
the limit given by the Vaughan-Preston gap dVaughan & Prestonl 
1980). The significant variability of active stars is produced 
by large magnetic features rotating with the star, that also in- 
duce large short-term RV variations. Indeed, all stars classi- 
fied as active exhibit a minimum short-term RV variability of 
10ms _I . Combined with imperfect sampling, this does not al- 
low us to properly estimate the RV - R' HK correlation slope with 
the method described here. 

There are also four inactive stars that seem to be out- 
liers: HD 85390, HD 192310, HD 72673 and HD 144628, which 
present no RV - R' HK correlation (small circles near 5200 K in 
Fig- HE). The first two st ars host multi-planet systems with long- 
period planets (see Mordasin i et al.l201 UlPepd201 lb . The long- 
period planets are not fully covered yet in phase and it is likely 
that the orbital solutions, and thus the RV residuals, will change 
when more points become available. Therefore, this can explain 
the observed lack of correlation. The third star presents a large 
peak in the periodogram at the rotational period of the star. We 
are therefore again confronted to short-term activity induced by 
magnetic features rotating with the star. In the case of an impor- 
tant signal related to short-term activity, even if we select only 3- 
month bins with more than 3 measurements, the sampling must 
be very dense to reduce the effect of activity. This is not the case 
for HD 72673, which can explain the unexpected value of the RV 

- R' HK correlation slope. 

Once active stars are removed from the plot (as well as the 
four stars mentioned above), we obtain a much cleaner relation 
between the RV - R^ correlation slope and effective tempera- 
ture, as shown in Fig. |T8] (bottom panel). The fitted linear model 
is given by: 



C 



RV 



(20.17 + 0.27) + (2.689 + 0.097) • 10" 2 ■ f efl 



(2.6 + 1.8)- 10" 6 -f 2 



eff 



(20.2 + 1.3) -[Fe/H] 



+ (26.3 + 2.5) ■ [Fe/H] 2 



(13) 



The different variables are the same as in Eq. [9] above. 
Comparing the two equations, we see that this second method, 
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Fig. 18. Top: Correlation slope between RV and R' nK , as a func- 
tion of effective temperature. Circles and squares are used to 
indicate inactive stars (max(log R' HK ) < -4.75) and active stars 
(max(log R' HK ) > -4.75), respectively (the maximum is calcu- 
lated on the 3 -month bins). In addition, three different sym- 
bol sizes are used to show the Pearson correlation coefficient 
R of the RV - R' nK correlation: small for R < 0.5, medium for 
0.5 < R < 0.75 and large for R > 0.75. Bottom: Same figure but 
with inactive stars only. 



which focuses on a smaller sample of carefully-selected stars, 
provides a T e ff dependence that is very similar to the previous 
one. Indeed, all the obtained polynomial coefficients are compat- 
ible with each other at the 2<x level. We conclude that both meth- 
ods confirm that 1) magnetic cycles do induce RV variations in 
solar-type stars, and 2) the sensitivity of RV measurements to a 
given magnetic cycle amplitude decreases towards cooler stars. 

5. Discussion and conclusion 

5.1. Calibration of R' HK taking into account metallicity 

In this paper we provide an updated calibration of the R' HK in- 
dex, which takes into account the effects of metallicity on stellar 
bolometric luminosity. Based on a calibration of effective tem- 
perature as a function of B - V color and metallicity, a correction 
factor to the usual R' UK definition can be applied and removes the 
observed trend of increasing R' values with decreasing metal- 



licity. The correction yields a sharper distribution of mean R' HK 
values for our sample stars, confirming that stars of different 
metallicities were exhibiting an artificial spread in R' HK with the 
uncorrected calibration. 



5.2. Statistics of magnetic cycles 

The present survey has allowed us to revisit some of the prop- 
erties of magnetic cycles in solar-type stars. To the best of our 
knowledge, most of the information previously available on this 
subject comes from the Mt Wilson project, and in particular the 
global analysis of lBaliunas et al.l (fl995) (see however Saar 2009, 
for a recent study). Two main difficulties arise when aiming at 
characterizing magnetic cycles: the required long-term preci- 
sion on Ca II H&K measurements, and the necessity to monitor 
many stars over decadal timescales. It turns out that large radial- 
velocity planet-search surveys, which started in the 1990s, can 
actually provide both, since they are based on high-resolution 
echelle spectrographs and regularly monitor hundreds of solar- 
type stars over many years. Several surveys have already yielded 
mean activity levels for many stars, but time series of activity 
measurements have been rarely discussed. 

The HARPS data presented here are of excellent qual- 
ity in terms of precision on R' HK , but they span "only" 7-8 
years, compared to 25 years for the published Mt Wilson data 
dBahunas et aljfl995l) . However, as shown in Sect. 13.41 this did 
not prevent us from detecting magnetic cycles with estimated 
periods up to ~15 yr, and long-term drifts of undetermined pe- 
riodicity. Obviously, the precision on fitted periods quickly de- 
creases when going beyond the time span of observations, but 
even uncertainties of 25-50% on the period still provide useful 
information on cycle properties. 

An important simplification adopted in this work is the 
assumption of perf ectly sinusoidal magnetic cycles (as in 
Baliu nas et alJll995h . We know from the Sun that they are in 
fact not really sinusoidal for any given cycle number, and more- 
over their amplitudes and shapes may vary significantly from 
one cycle to the next. However, our survey usually covers only 
one or two cycles, so that we are not able to study the evolu- 
tion of the magnetic cycle in any given star. We rather obtain a 
random snapshot of cycles in many stars at the same time. With 
these limitations in mind, we can summarize the main results of 
this work as follows: 

- The global occurrence of magnetic cycles with semi- 
amplitude A > 0.04 is 61±8% among older solar-type stars. 
Correspondingly, 39+8% of older solar-type stars show a flat 
activity record. These numbers are roughly compatible with 
the results of Baliunas et al. (1995), although these authors 
do not explicitly provide comparable estimates and use a 
different classification scheme. The picture is markedly dif- 
ferent for subgiants, for which the occurrence of magnetic 
cycles drops to only 12+5%. 

- The distribution of magnetic cycle periods shows a broad 
maximum between 2-10 yr, and a sharp decrease towards 
longer periods. This may be due in part to the limited du- 
ration of the survey, but the number of observed long-term 
trends does not seem to be sufficient to produce a large num- 
ber of long-period cycles. The decrease towards longer pe- 
riods therefore seems to be real. There may be some dis- 
crepancies between our period distribution and the one of 
Baliunas et al. ( 1995), in the sense that we detect more short- 
period cycles than they do. It must be noted, however, that 
short-period cycles are often less well-defined (more disper- 
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sion and evolution with time) than longer-period ones, and 
the discrepancy may be explained by a different classifica- 
tion s cheme combin e d w ith a longer duration of observa- 
tions (Baliunas et al.l (1 19951) classify as "variable" the stars 
with poorly defined periodicities). On the other hand, both 
studies seem to agree on the low number of cycles with P > 
15 yr. This confirms that the limited duration of our survey 
does not significantly affect the derived overall properties of 
magnetic cycles. 

The distribution of magnetic cycle semi-amplitudes shows a 
broad peak between zero and 0.2 in R' HK - Given the obser- 
vational bias against low amplitudes, it seems that magnetic 
cycles significantly less pronounced than the solar one are 
the rule rather than the exception. 

We find no clear correlations between cycle periods and am- 
plitudes, nor between cycle periods and mean activity levels. 
The sometimes large error bars on cycle periods prevent us 
from going into further details. The occurrence and stability 
in time of short-period cycles (P < 5 yr) should be further in- 
vestigated, together with methods to disentangle them from 
evolution of active regions on rotational timescales. 
We find a clear trend of increasing cycle amplitudes with 
increasing mean activity levels. In particular, no stars with 
mean logR' HK a -5.0 show cycles as pronounced as the Sun. 
We also find that the upper envelope of cycle amplitudes 
shows a peak in the temperature range 5000-5600 K, while 
hotter stars do not exhibit very large amplitudes. 
The distribution of mean log R' HK values is markedly differ- 
ent between G and K dwarfs, in the sense that G dwarfs tend 
to accumulate around -5.0, while K dwarfs are more spread 
between -4.7 and -5.0. A possible explanation is the slower 
decrease in activity with age in K dwarfs, which need more 
than ~6-8 Gyr to reach a mean level of -5.0. 
Non-cycling stars are only found at mean activity levels 
lower than about -4.95, and almost only among G dwarfs. 
Cycling stars are preferentially found at higher activity lev- 
els, although some of them are also found around -5.0. The 
number of well-observed K dwarfs below -4.95 is too low 
to properly study the occurrence of non-cycling stars in this 
region of parameter space. 

A possible interpretation of these results is that magnetic cy- 
cle amplitudes decrease with age, and tend to disappear when 
stars reach a mean activity level of about -5.0. This is also 
supported by the quasi-absence of cycles in subgiants. 
The Sun, with a mean logR' UK of -4.91, is still in the ac- 
tivity range where most stars exhibit magnetic cycles. One 
can conjecture that solar cycles will slowly disappear as the 
mean activity decreases to -5.0, possibly within a few Gyr. 
It is difficult to estimate from our data how frequent peri- 
ods of extended minima like the Maunder minimum are for 
any given star, since we only have instantaneous snapshots 
of the state of many stars. What can be said is that the activ- 
ity level of the Sun was almost certainly below -4.95 during 
the Maunder minimum since we do not find any non-cycling 
stars above that threshold. Considering the well-established 
decrease in activity with age in solar-type stars, it is more 
likely that the "clump" of G dwarfs around -5.0 is mainly 
the signature of an age effect, and not a population of stars in 
temporary Maunder-minimum states. However, an unknown 
proportion of these stars are probably in such a state. This 
also supports the view that Maunder-minimum stars can- 
not be found on the basis of their mean activity level only 
(Wright 2004), and should rather be identified by the demon- 
strated absence of a magnetic cycle. 
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Fig. 19. Effective impact of magnetic cycles on RV for all sample 
stars with a detected cycle, shown as a function of T s s . 



5.3. Effects of magnetic cycles on spectral lines and radial 
velocities 

It has long been hypothesized that magnetic cycles have an 
impact on photpspheric spectral lines (e. g. iDravinsI [1985; 
iGrav et al.|[l992l:lLindegren & Dravinsll2003l) . Here we provide 
a global assessment of magnetically-induced variability in four 
quantities derived from the HARPS cross-correlation function 
(RV, FWHM, contrast and BIS), which is constructed as a "mas- 
ter" spectral line built from thousands of mainly weak to moder- 
ate individual Fe lines. 

We observe long-term correlations between R' HK , RV, 
FWHM, contrast and BIS in many stars of our sample, with 
varying sensitivities between stars of different spectral types and 
various behaviors for different quantities. Globally speaking, for 
a given variation in R' HK , we find that RV and BIS are more sen- 
sitive in G dwarfs than in K dwarfs, whereas the opposite is true 
for the FWHM. The contrast sensitivity is roughly constant with 
spectral type, but varies significantly with metallicity. Overall, 
sensitivity variations with temperature are more important than 
with metallicity for RV, FWHM and BIS. The sensitivity of RV 
to magnetic cycles is almost zero at Tes ~ 4800 K, while the 
same is true for FWHM at T e ff w 6200 K. However, the effective 
impact of a magnetic cycle for any given star is obtained by mul- 
tiplying these sensitivities with the actual R' HK semi-amplitude, 
which may be large and therefore have a non-negligible effect 
even if the sensitivity is low. We show in Fig. [19] the actual im- 
pact of magnetic cycles on the RV data for all stars in the sample 
with a detected cycle. For each star the RV sensitivity is com- 
puted from the model given in Eq.|9j and the result is multiplied 
by the measured R' HK semi-amplitude A. We see that induced 
RV semi-amplitudes can reach up to 11 m s in the worst case. 
More typical values are between and 3 ms _1 (there are unaf- 
fected stars at all spectral types). Interestingly, the model pre- 
dicts an anti-correlation between R' nK and RV (i.e. a negative 
RV semi-amplitude in Fig. [19]) for stars cooler than ~4800 K. At 
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least in one case this surprising result se ems to be c onfirmed by 
an in-depth analysis: the star HD 85512 (Pepe 201 1). Looking at 
Fig -HI we also see that despite the drop in sensitivity towards 
cool stars, one can still have RV effects of several m s because 
of the occurrence of strong magnetic cycles in K dwarfs. 

The previous study that is most directly comparabl e to the 
results presented here is the one bv lSantos et al.1 (I2010I) . These 
authors obtained long-term HARPS data on eight stars known 
to have clear magnetic cycles, and studied the behavior of CCF 
parameters in the same way as we do here. Clear correlations 
were found between the S index, FWHM and BIS data, while 
the contrast was anti-correlated to the S index. These findings 
are in excellent agreement with the present results. On the other 
hand, attempts to find correlations with the RV data remained 
inconclusive, or actually RV effects were too small to be de- 
tected. The sensitivity model derived in this paper provides a 
natural explanation for this lack of cor relations: it turns ou t that 
the five cycling stars actually studied by S an tos et all ([20 1 ) span 
a very narrow range of effective temperatures (200 K), centered 
at 5000 K. In this regime RV sensitivities are close to zero, and 
the S index variations were not large enough to induce a measur- 
able RV effect. The only hot star in the sample, HD 216385, is 
a magnetically-constant star that was chosen as a standard. Two 
other, higher-temperature stars in the sample could not be stud- 
ied because of strong short-term variability and binarity. We con- 
clude that it is only bad luck that prevented Sa ntos et al.l ([2010) 
from finding correlations between RV and S index. 

In general, the RV sensitivities to magnetic cycles are such 
that the induced RV effects can be of the same order of magni- 
tude and period as the signal of Jupiter-like planets orbiting at 
several AUs from the star. Therefore, it is critical to check the 
behavior of the chromospheric activity index before announc- 
ing any long-period planet with a RV semi-amplitude lower than 
~15 mf 1 . It would be extremely useful for all planet-search 
teams to measure chromospheric activity (even at lower preci- 
sion) in all stars exhibiting such low-amplitude, long-period RV 
signals, and to check whether any correlations are found with the 
magnetic cycle. This obviously also applies to already published 
planets, e.g. the Ju piter-analog candidate around HD 154345 
(Wright et al. 2008). If simultaneous activity measurements are 
available, then the present work shows that they can be easily 
used to correct RV data from magnetic cycle effects. Two ap- 
proaches can be used: either one includes an additional sinu- 
soidal component in the Keplerian model, fixing its period and 
phase to the ones of the magnetic cycle, or one can use the sensi- 
tivity model presented here to a priori correct the RV data based 
on a fit of the magnetic cycle on the activity data. Several exam- 
ples of these correction procedures and their accuracies can be 
found in Dumusque et al. (2 011al) . 

The sensitivity model also allows us to estimate the RV ef- 
fect that is expected in the Sun due to the solar magnetic cycle. 
With r eff = 5770 K, [Fe/H] = and A = 0.26, we obtain a RV 
semi-amplitude of 8.2 ms -1 and thus a peak-to-peak effect of 
about 16 ms" 1 . This i s in qualitative agreement w ith both the- 
oretical estimates (e.g.[Lindegren & Dravin sll2003l) and a previ - 
ous attempt to actually measure it (iDeming & Plvmata fl994). 
We note that the effect is likely to depend on the particular set of 
spectral lines being considered, so that the comparison with lit- 
erature values is not straightforwa rd. This coul d potentially ex- 
plain the null result of McMillan et alj d!993l) . A more recent 
attempt to predict this same effect and its impact on exoplanet 
searches was published by Meunier et al.l ((2010). They obtain 
a full amplitude over the solar cycle of about 8 ms -1 , which 
is lower than our result but can also be considered as qualita- 



tively in agreement given the differences in the details between 
our observations and their simulations. We note in passing that 
Meu nier et al.l (2010) went on to simulate planet detection limits 
in the presence of stellar noise, and concluded that the effects 
of magnetic cycles would make the detection of Earth-like plan- 
ets impossible. We emphasize here that these long-term effects 
are actually rather easy to model and correct using simultaneous 
R' HK measurements. Their planet detection limits are therefore 
likely far too pessimistic. 

Finally, we briefly discuss theoretical explanations for the 
observed effects of magnetic cycles on spectral lines. As al- 
ready argued by several authors, the main phenomenon at play 
here is likely to be the "freezing" effect of magnetic fields on 
gran ular motions in the outer convective zone of solar-type stars 
(e.g. lLivingstonll982tlBrandt & Solankil 19901) . The stronger the 
magnetic field, the slower the convective motions. As a conse- 
quence, the usual convective blueshift of spectral lines , caused 
by the temperature-velocity c orrelation in granules dPravinsi 
11982b iKaisig & Durrantl Il982l) . is decreased in active regions 
compared to the quiet photosphere. Therefore, stellar lines will 
appear redshifted when the average (or local) magnetic field 
strength is higher. This may explain the observed RV - R' HK pos- 
itive correlation over the magnetic cycle. However, other effects 
may also play a role, suc h as changes in l arge-scale surface flows 
over the magnetic cycle (Makarov 2010), or perhaps changes in 
global mean effective temperature as the total spot/plage cover- 
age varies (sensitivity of spectral lines to temperature). 

Differential effects as a function of effective temperature or 
metallicity are more difficult to explain. Convective motions are 
slower in cooler stars, diminishing the amount of convective 
blueshift and therefore also reducing possible variations with 
magnetic field strength. The strong increase in FWHM sensi- 
tivity towards cooler temperatures should also be considered to 
help us understand the effect. The direct impact of magnetic 
fields on spectral lines, through Zeeman broadening, should be 
investigated. 3D hydrodynamical simulations of stellar atmo- 
spheres including magnetic fields are probably necessary to bet- 
ter understand all these observational findings. On the observa- 
tional side, the lists of spectral lines that are used to build the 
cross-correlation function for precise RV measurements should 
be carefully checked and optimized to either minimize undesir- 
able RV effects, or actually maximize them in order to better un- 
derstand the underlying physics affecting line shapes over mag- 
netic cycles. 

In conclusion, the present study has allowed us to make im- 
portant progress in the knowledge of magnetic cycles in solar- 
type stars in general, and in their effects on spectral lines in par- 
ticular. Not only are these results hopefully useful in their own 
right for stellar physics, but they are also of great help for the 
search for exoplanets with the RV method. Understanding stel- 
lar RV "jitter" is the single most important challenge in the quest 
for Earth-like planets in the habitable zones of solar-type stars. 
The excellent quality of the HARPS data, and the different di- 
agnostics available from the spectra, make it possible to monitor 
stars with unprecedented scrutiny. Further studies of the impact 
of magnetic activity at all timescales are ongoing and will al- 
low us to disentangle even better stellar activity and planetary 
signals. 
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0.0177 


8607 40 


1210 


-4.722 


-4.713 


0.1623 


0.1069 


8978 A 61 


2340 


-4.875 


-4.876 


0.0459 


0.0364 


9830 205 


2747 


-4.986 


-4.996 


0.0565 


0.0217 


1479 16 


2355 


-5.015 


-5.016 


0.0382 


0.0149 


1835 72 


2698 


-4.898 


-4.907 


0.1052 


0.0588 


2673 65 


2464 


-4.968 


-4.961 


0.0668 


0.0233 


3121 20 


2709 


-5.048 


-5.048 


0.0219 


0.0121 


3524 60 


2694 


-4.981 


-4.981 


0.0385 


0.0197 


8747 43 


2554 


-4.991 


-4.992 


0.0167 


0.0121 


2342 31 


3009 


-5.012 


-5.010 


0.0753 


0.0333 


2516 56 


2320 


-4.954 


-4.963 


0.1222 


0.0537 


5390 65 


2689 


-4.969 


-4.967 


0.0631 


0.0401 


5512 175 


2705 


-4.936 


-4.937 


0.1143 


0.0317 


6140 17 


2690 


-4.838 


-4.844 


0.1002 


0.0506 


6065 24 


3008 


-4.872 


-4.877 


0.1245 


0.0621 


9454 48 


1468 


-4.685 


-4.688 


0.1353 


0.0906 


0156 79 


2665 


-4.973 


-4.973 


0.0131 


0.0111 


0711 18 


2306 


-4.974 


-4.972 


0.0837 


0.0269 


0812 22 


2307 


-4.981 


-4.994 


0.1170 


0.0286 


3083 66 


2279 


-4.941 


-4.933 


0.1627 


0.0443 


4151 20 


2703 


-4.960 


-4.981 


0.0915 


0.0233 


6700 152 


3011 


-4.967 


-4.968 


0.0199 


0.0178 


8281 53 


2227 


-4.914 


-4.911 


0.0670 


0.0487 


00508 32 


2284 


-5.001 


-5.008 


0.0776 


0.0297 


01930 33 


2261 


-5.016 


-5.019 


0.0497 


0.0304 


06116 103 


2658 


-5.005 


-5.007 


0.0243 


0.0143 


06275 17 


2649 


-4.870 


-4.891 


0.1956 


0.0709 


08309 20 


2647 


-5.003 


-5.000 


0.0341 


0.0163 


09200 118 


2866 


-4.976 


-4.981 


0.0771 


0.0335 


14613 96 


1768 


-5.059 


-5.063 


0.0545 


0.0371 



■4.993 


0.0351 


29.0 ± 3.3 


5.12 ±0.59 


■4.953 


0.0099 


16.7 ± 2.6 


2.92 ± 0.34 


■4.820 


0.1804 


21.3 + 3.2 


3.07 ± 0.35 


■4.916 


0.2620 


41.4 + 4.7 


6.38 ± 0.74 


■5.024 


0.0383 


38.8 + 3.8 


7.05 ± 0.81 


■4.883 


0.1998 


39.8 + 4.7 


5.85 ± 0.67 


■4.999 


0.0216 


24.1 +3.0 


4.28 ± 0.49 


■4.766 


0.5957 


29.8 ± 4.3 


3.90 ± 0.45 


■4.759 


0.1582 


19.9 + 3.3 


2.63 ± 0.30 


■4.883 


0.2688 


40.2 + 4.8 


5.87 ± 0.68 


■4.902 


0.1185 


41.8 + 4.9 


6.20 ± 0.71 


■4.934 


0.1711 


37.1 +4.1 


6.08 ± 0.70 


■4.905 


0.0495 


17.7 + 2.8 


2.87 ± 0.33 


■4.846 


0.1438 


22.3 + 3.2 


3.32 ±0.38 


■5.025 


0.0065 


35.4 + 3.6 


6.45 ± 0.74 


■5.080 


0.0848 


49.5 + 4.4 


9.19 ± 1.06 


■4.920 


0.0474 


21.6 + 3.0 


3.52 ± 0.41 


■4.733 


0.2037 


34.5 + 6.8 


3.32 ±0.38 


■4.935 


0.2094 


36.3 + 4.0 


5.98 ± 0.69 


■4.958 


0.1585 


42.2 + 4.4 


6.95 ± 0.80 


■4.749 


0.0903 


19.5 + 3.3 


2.54 ± 0.29 


■4.932 


0.2842 


41.3 + 4.5 


6.57 ± 0.76 


■5.032 


0.0424 


41.4 + 4.0 


7.54 ± 0.87 


■4.945 


0.0208 


23.6 + 3.1 


3.95 ± 0.45 


■4.957 


0.1894 


47.2 + 5.3 


7.08 ± 0.81 


■4.913 


0.0979 


21.4 + 3.0 


3.45 ± 0.40 


■4.970 


0.0717 


34.0 + 3.7 


5.86 ± 0.67 


■4.981 


0.0317 


39.6 + 4.1 


6.84 ± 0.79 


■4.981 


0.0773 


27.4 + 3.2 


4.76 ± 0.55 


■5.000 


0.0361 


31.8 + 3.5 


5.66 ± 0.65 


■4.772 


0.2275 


26.2 + 3.8 


3.52 ± 0.41 


■4.923 


0.0779 


14.9 + 2.6 


2.57 ± 0.30 


■5.037 


0.0949 


54.9 + 6.2 


8.05 ± 0.93 


■4.881 


0.1910 


17.1 +2.8 


2.69 ±0.31 


■4.725 


0.1495 


29.9 + 4.7 


3.54 ± 0.41 


■4.867 


0.0566 


18.8 + 2.9 


2.88 ± 0.33 


■4.944 


0.1439 


35.8 + 3.9 


5.98 ± 0.69 


■4.997 


0.0798 


26.6 + 3.2 


4.72 ± 0.54 


■4.923 


0.1803 


36.5+4.1 


5.91 ± 0.68 


■4.975 


0.0910 


40.2 + 4.1 


6.86 ± 0.79 


■5.057 


0.0343 


19.2 + 2.6 


3.72 ± 0.43 


■4.981 


0.0464 


19.8 + 2.8 


3.49 ± 0.40 


■4.994 


0.0188 


16.5 + 2.6 


3.06 ± 0.35 


■4.999 


0.1184 


50.9 + 5.6 


7.68 ± 0.88 


■4.802 


0.5874 


36.7 + 5.3 


4.59 ± 0.53 


■4.973 


0.0790 


44.9 + 4.7 


7.37 ± 0.85 


■4.956 


0.2035 


50.9 + 7.0 


6.41 ± 0.74 


■4.831 


0.1836 


39.0 + 5.4 


4.98 ± 0.57 


■4.856 


0.1275 


40.7 + 5.4 


5.37 ± 0.62 


■4.695 


0.1606 


21.1 +3.7 


2.51 ±0.29 


■4.973 


0.0092 


27.0 + 3.2 


4.67 ± 0.54 


■4.984 


0.1180 


43.1 +4.4 


7.33 ± 0.84 


■4.831 


0.5507 


34.1 +4.4 


4.85 ± 0.56 


■4.942 


0.2273 


46.2 + 5.3 


6.81 ±0.78 


■4.951 


0.1561 


33.1 +3.7 


5.58 ± 0.64 


■4.965 


0.0141 


20.6 + 2.9 


3.55 ± 0.41 


■4.938 


0.1120 


33.5 + 3.8 


5.57 ± 0.64 


■4.903 


0.3236 


39.5 + 4.5 


6.09 ± 0.70 


■4.906 


0.3137 


43.0 + 5.1 


6.25 ± 0.72 


■5.002 


0.0335 


33.1 +3.5 


5.90 ± 0.68 


■4.893 


0.3246 


41.6 + 5.0 


6.05 ± 0.70 


■5.014 


0.0547 


29.9 + 3.3 


5.39 ± 0.62 


■4.950 


0.2118 


43.0 + 4.6 


6.92 ± 0.80 


■5.053 


0.0553 


34.1 +3.5 


6.38 ± 0.73 
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Table 2. Continued. 



C. Lovis et al.: The HARPS search for southern extra-solar planets 



Name 



Span 
[days] 



log^ K 
mean 



log^ K 
median 



raw 



short-term 



'cycle 
[days] 

1820!"' 

1548 -|) 
1947 + |* 

6000!~ 

1146 + ^ 2 

1044 + « 
1419+336 

2097 +o ° 

**" ' -691 

1041 + 5|i 

5706!™,! 

2180+5" 

2770+^' 

4103+ 1321 

2803! 2 *f 

2416!|« 4 

8000!" 

2702!™ 

8000!" 

1613!™ 

8000!™, , 

1221+ 481 

2016+ 355 

-262 

999 1 +>w 
-431 

3792!™° 

764 +2 f 

1076+* 3 

2080!^ 2 

987!" 

21 32 -232 

1520!"' 

1719+217 

713™ 
2168™* 

1667 + ^< 
1336!^° 



y 



'rot 

[days] 



Age 
[Gyr] 



14747 


38 


1473 


-4.829 


-4.859 


0.2596 


0.0767 


15617 


129 


2639 


-4.992 


-4.995 


0.0241 


0.0198 


15674 


39 


2252 


-4.919 


-4.924 


0.0484 


0.0295 


19638 


24 


3000 


-4.934 


-4.932 


0.0417 


0.0196 


25072 


24 


2174 


-4.941 


-4.952 


0.1343 


0.0555 


25881 


26 


3003 


-4.872 


-4.867 


0.0795 


0.0289 


29642 


54 


1510 


-4.971 


-4.973 


0.0488 


0.0363 


32648 


28 


2624 


-4.874 


-4.890 


0.1187 


0.0528 


34664 


28 


1066 


-4.871 


-4.858 


0.0823 


0.0263 


36352 


123 


2543 


-4.986 


-4.986 


0.0107 


0.0098 


36713 


41 


2203 


-4.787 


-4.771 


0.1429 


0.0710 


38549 


21 


2622 


-4.834 


-4.819 


0.1370 


0.0639 


37388 


64 


2118 


-4.879 


-4.867 


0.1461 


0.0723 


44628 


53 


2105 


-5.004 


-5.003 


0.0508 


0.0280 


46233 


38 


2157 


-4.923 


-4.923 


0.0793 


0.0287 


48303 


31 


2164 


-4.667 


-4.657 


0.2140 


0.1017 


54088 


107 


1866 


-5.034 


-5.037 


0.0267 


0.0174 


54577 


118 


2550 


-4.974 


-4.981 


0.0750 


0.0289 


57172 


81 


2099 


-4.990 


-5.001 


0.0798 


0.0228 


57830 


55 


2526 


-4.822 


-4.824 


0.1164 


0.0655 


61098 


65 


1866 


-4.946 


-4.952 


0.0500 


0.0282 


72513 


41 


1081 


-4.777 


-4.777 


0.0751 


0.0503 


77565 


26 


1684 


-4.927 


-4.949 


0.1218 


0.0248 


88559 


35 


1153 


-4.811 


-4.826 


0.1981 


0.0952 


92310 


152 


2704 


-4.996 


-5.009 


0.1060 


0.0258 


97823 


19 


1056 


-4.719 


-4.751 


0.2345 


0.0377 


99960 


30 


1864 


-4.982 


-4.978 


0.0488 


0.0174 


03432 


36 


1405 


-4.816 


-4.838 


0.1904 


0.0641 


04313 


66 


1548 


-4.993 


-4.997 


0.0323 


0.0233 


04941 


34 


2180 


-4.962 


-4.971 


0.0695 


0.0238 


07129 


79 


1876 


-4.898 


-4.907 


0.0890 


0.0400 


09100 


46 


1916 


-4.806 


-4.806 


0.1025 


0.0641 


15152 


167 


2740 


-4.871 


-4.878 


0.0837 


0.0516 


22237 


18 


1115 


-5.015 


-5.025 


0.0575 


0.0246 


22595 


26 


1053 


-4.806 


-4.833 


0.2024 


0.0833 


22669 


47 


1403 


-4.858 


-4.860 


0.0715 


0.0526 



-4.794 
-4.993 
-4.902 
-4.932 
-4.859 
-4.867 
-4.968 
-4.845 
-4.871 
-4.986 
-4.921 
-4.851 
-4.877 
-5.017 
-4.918 
-4.674 
-4.977 
-4.959 
-4.840 
-4.822 
-4.801 
-4.799 
-4.895 
-4.811 
-4.931 
-4.645 
-4.992 
-4.838 
-5.000 
-4.957 
-4.899 
-4.788 
-4.867 
-4.993 
-4.796 
-4.858 



0.3770 
0.0200 
0.0852 
0.0484 
0.2797 
0.1241 
0.0464 
0.2082 
0.1588 
0.0060 
0.5547 
0.1759 
0.1844 
0.0797 
0.1476 
0.3000 
0.1550 
0.1294 
0.4855 
0.1487 
0.4861 
0.1205 
0.1710 
0.3193 
0.2049 
0.4882 
0.0857 
0.2825 
0.0371 
0.0892 
0.1069 
0.1486 
0.0876 
0.0918 
0.2280 
0.0654 



36.0 + 5.2 
33.9 + 3.6 

25.5 + 3.3 
12.0 ± 2.3 
42.0 ± 5.9 

15.7 + 2.7 

48.0 ± 5.3 

28.1 +3.7 

23.2 + 3.2 

25.0 + 3.1 
45.4 + 5.6 

28.4 + 3.7 

40.6 + 5.0 

48.2 + 4.7 

23.8 + 3.2 

28.9 + 5.6 

42.6 + 4.4 
45.8 + 5.0 

32.7 + 4.1 

24.3 + 3.5 

22.5 + 3.4 

28.3 + 3.9 

28.4 + 3.5 

39.1 +6.2 

43.7 + 4.9 

21.0 + 4.0 

25.2 + 3.1 

29.5 + 3.8 

33.1 +3.5 

45.3 + 4.9 

17.6 + 2.8 

37.6 + 6.2 

41.8 + 5.6 

50.7 + 5.7 
24.1 +3.5 
17.5 + 2.9 



4.49: 
5.99: 
4.08: 
2.35: 
5.23: 
2.45: 
7.27: 
4.18: 
3.58: 
4.37: 
6.36: 
4.26: 
5.79: 
8.21 : 
3.87: 
2.82: 
7.20: 
7.18: 
4.78: 
3.51 : 
3.16: 
3.95: 
4.50: 
4.45: 
6.69: 
2.27: 
4.43: 
4.35: 
5.88: 
7.15: 
2.83: 
4.12: 
5.49: 
7.54: 
3.36: 
2.66: 



0.52 
0.69 
0.47 
0.27 
0.60 
0.28 
0.84 
0.48 
0.41 
0.50 
0.73 
0.49 
0.67 
0.95 
0.45 
0.33 
0.83 
0.83 
0.55 
0.40 
0.36 
0.45 
0.52 
0.51 
0.77 
0.26 
0.51 
0.50 
0.68 
0.82 
0.33 
0.47 
0.63 
0.87 
0.39 
0.31 
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C. Lovis et al.: The HARPS search for southern extra-solar planets 
Table 3. Results for stars with no significant magnetic cycle. 



Name 


N 

1 v meas 


Span 


l°g*HK 


log^ K 


^HK> 


^hk) 


Alim 


Prob 


*rot 


Age 






[days] 


mean 


median 


raw 


short-term 




(A < 0.04) 


[days] 


[Gyr] 


HD224789 


32 


1328 


-4.436 


-4.438 


0.1638 


0.1460 


0.5968 


0.000 


11.8 ±4.7 


0.68 ± 0.08 


HD55 


10 


1311 


-4.893 


-4.891 


0.0965 


0.0283 


0.6010 


0.000 


45.2 ± 6.4 


5.59 ± 0.64 


HD361 


16 


2134 


-4.712 


-4.702 


0.1259 


0.0878 


0.5226 


0.000 


14.1 ±2.9 


1.73 ±0.20 


HD967 


38 


2556 


-4.987 


-4.986 


0.0218 


0.0192 


0.0468 


0.833 


26.3 ± 3.2 


4.60 ± 0.53 


HD1320 


6 


2510 


-4.904 


-4.879 


0.0710 


0.0201 


0.6010 


0.000 


23.4 ± 3.2 


3.74 ± 0.43 


HD1388 


57 


2647 


-4.979 


-4.979 


0.0160 


0.0146 


0.0517 


0.917 


19.7 ±2.8 


3.47 ± 0.40 


HD3569 


7 


2081 


-4.938 


-4.965 


0.1808 


0.0430 


0.6010 


0.000 


42.7 ± 4.7 


6.75 ± 0.78 


HD3823 


33 


2266 


-5.017 


-5.015 


0.0135 


0.0118 


0.0375 


1.000 


15.8 ±2.5 


3.08 ± 0.35 


HD4308 


63 


2178 


-4.999 


-4.997 


0.0230 


0.0190 


0.0364 


1.000 


28.0 ±3.2 


4.96 ± 0.57 


HD4307 


31 


1856 


-5.069 


-5.068 


0.0147 


0.0139 


0.0443 


0.917 


22.2 ± 2.8 


4.26 ± 0.49 


HD6348 


8 


2559 


-4.983 


-4.990 


0.0389 


0.0236 


0.5749 


0.000 


42.3 ± 4.3 


7.22 ± 0.83 


HD6673 


9 


1142 


-4.818 


-4.815 


0.0804 


0.0292 


0.6010 


0.000 


37.7 ±5.3 


4.82 ± 0.56 


HD6735 


9 


2560 


-4.890 


-4.891 


0.0260 


0.0134 


0.3863 


0.000 


14.2 ± 2.6 


2.35 ± 0.27 


HD7134 


16 


1767 


-4.980 


-4.980 


0.0133 


0.0086 


0.0637 


0.750 


18.6 ±2.7 


3.29 ± 0.38 


HD7449 


82 


2633 


-4.860 


-4.861 


0.0477 


0.0407 


0.0642 


0.667 


13.5 ±2.6 


2.16 ±0.25 


HD8638 


35 


1827 


-4.996 


-4.996 


0.0216 


0.0165 


0.0726 


0.750 


30.4 ± 3.4 


5.39 ± 0.62 


HD9796 


12 


2112 


-4.920 


-4.919 


0.0658 


0.0239 


0.4264 


0.000 


37.1 ±4.1 


5.98 ± 0.69 


HD9782 


27 


2217 


-4.964 


-4.963 


0.0179 


0.0173 


0.1053 


0.000 


18.1 ±2.7 


3.16 ±0.36 


HD 10002 


8 


2088 


-5.051 


-5.049 


0.0230 


0.0117 


0.6010 


0.000 


48.5 ± 4.5 


8.72 ± 1.00 


HD 10700 


153 


2650 


-5.017 


-5.017 


0.0089 


0.0084 


0.0086 


1.000 


37.4 ±3.8 


6.75 ± 0.78 


HD11226 


28 


1766 


-4.998 


-4.998 


0.0114 


0.0108 


0.0597 


0.750 


16.6 ± 2.6 


3.09 ± 0.36 


HD11505 


12 


1392 


-5.024 


-5.026 


0.0136 


0.0123 


0.1396 


0.000 


26.3 ±3.1 


4.78 ± 0.55 


HD11964A 


60 


2514 


-5.154 


-5.154 


0.0152 


0.0135 


0.0184 


1.000 


53.9 ± 4.4 


10.69 ± 1.23 


HD12387 


14 


1092 


-5.017 


-5.004 


0.0544 


0.0364 


0.5131 


0.000 


27.4 ± 3.2 


4.94 ± 0.57 


HD 12345 


14 


1729 


-5.017 


-5.013 


0.0289 


0.0168 


0.1242 


0.167 


39.5 ± 3.9 


7.12 ±0.82 


HD12617 


12 


1544 


-4.694 


-4.717 


0.2616 


0.1161 


0.6010 


0.000 


30.7 ± 5.9 


3.00 ± 0.35 


HD 14374 


17 


1805 


-4.660 


-4.671 


0.1760 


0.1314 


0.6010 


0.000 


20.6 ± 3.8 


2.31 ±0.27 


HD 14747 


12 


2092 


-5.004 


-5.001 


0.0275 


0.0246 


0.2822 


0.000 


30.8 ± 3.4 


5.48 ± 0.63 


HD 16297 


10 


2052 


-4.705 


-4.710 


0.1561 


0.0754 


0.6000 


0.000 


24.7 ±4.1 


2.96 ± 0.34 


HD16417 


76 


2623 


-5.029 


-5.031 


0.0307 


0.0254 


0.0458 


0.833 


27.8 ± 3.2 


5.08 ± 0.58 


HD16714 


26 


2089 


-4.988 


-4.987 


0.0188 


0.0159 


0.0737 


0.417 


33.7 ±3.6 


5.92 ± 0.68 


HD17970 


8 


1739 


-5.071 


-5.065 


0.0225 


0.0141 


0.5599 


0.000 


49.7 ± 4.5 


9.10 ± 1.05 


HD 19034 


8 


2136 


-4.957 


-4.951 


0.0327 


0.0172 


0.5048 


0.000 


28.9 ± 3.4 


4.90 ± 0.56 


HD 19467 


15 


2170 


-5.015 


-5.015 


0.0149 


0.0122 


0.0697 


0.583 


27.3 ± 3.2 


4.92 ± 0.57 


HD20407 


8 


2524 


-4.943 


-4.942 


0.0189 


0.0081 


0.4619 


0.000 


17.6 ±2.7 


3.00 ± 0.35 


HD21209A 


9 


2707 


-4.876 


-4.883 


0.0872 


0.0322 


0.5814 


0.000 


43.9 ± 6.3 


5.36 ± 0.62 


HD21019 


13 


1800 


-5.122 


-5.122 


0.0166 


0.0151 


0.0735 


0.333 


38.4 ±3.5 


7.66 ± 0.88 


HD21938 


12 


2702 


-4.990 


-4.986 


0.0241 


0.0193 


0.2066 


0.000 


18.8 ±2.7 


3.37 ± 0.39 


HD22879 


54 


2687 


-4.996 


-4.997 


0.0140 


0.0119 


0.0259 


1.000 


14.2 ± 2.4 


2.83 ± 0.33 


HD23456 


20 


2201 


-4.939 


-4.937 


0.0216 


0.0154 


0.0930 


0.417 


9.4 ±2.1 


2.85 ± 0.33 


HD23249 


54 


1877 


-5.181 


-5.181 


0.0061 


0.0057 


0.0117 


1.000 


61.4 ±5.2 


11.56 ± 1.33 


HD23356 


13 


2228 


-4.774 


-4.769 


0.0748 


0.0447 


0.3493 


0.000 


34.8 ± 5.3 


4.19 ±0.48 


HD24331 


11 


2702 


-4.820 


-4.808 


0.1848 


0.0817 


0.6010 


0.000 


37.5 ±5.1 


4.89 ± 0.56 


HD24892 


13 


2749 


-5.124 


-5.123 


0.0146 


0.0111 


0.1421 


0.250 


45.0 ± 3.9 


8.95 ± 1.03 


HD25673 


13 


2123 


-4.971 


-4.995 


0.1436 


0.0562 


0.5434 


0.000 


43.0 ± 4.4 


7.17 ±0.83 


HD28471 


13 


1194 


-5.010 


-5.011 


0.0432 


0.0268 


0.3656 


0.000 


27.1 ±3.2 


4.87 ± 0.56 


HD28701 


19 


1330 


-5.021 


-5.017 


0.0230 


0.0123 


0.1494 


0.083 


27.5 ± 3.2 


4.98 ± 0.57 


HD28821 


12 


2706 


-4.984 


-4.985 


0.0314 


0.0140 


0.4125 


0.000 


29.9 ± 3.4 


5.24 ± 0.60 


HD30306 


15 


1454 


-5.052 


-5.053 


0.0262 


0.0199 


0.1436 


0.250 


41.4 ±3.9 


7.71 ±0.89 


HD30278 


21 


1514 


-5.024 


-5.021 


0.0235 


0.0165 


0.0892 


0.667 


39.9 ± 3.9 


7.23 ± 0.83 


HIP22059 


10 


1357 


-4.515 


-4.508 


0.0909 


0.0749 


0.6010 


0.000 


17.9 ±5.4 


1.25 ±0.14 


HD34449 


12 


2664 


-4.886 


-4.888 


0.0257 


0.0219 


0.3564 


0.000 


18.3 ±2.9 


2.88 ± 0.33 


HD31527 


166 


2720 


-4.972 


-4.972 


0.0127 


0.0122 


0.0140 


1.000 


20.8 ± 2.9 


3.61 ± 0.42 


HD31822 


45 


2353 


-4.885 


-4.886 


0.0216 


0.0206 


0.0415 


0.917 


15.1 ±2.6 


2.44 ± 0.28 


HD32724 


13 


2704 


-5.045 


-5.044 


0.0228 


0.0212 


0.2045 


0.000 


24.2 ± 2.9 


4.51 ±0.52 


HD33725 


10 


2720 


-4.970 


-5.003 


0.1040 


0.0107 


0.5416 


0.167 


41.5 ±4.3 


7.00 ± 0.81 


HD35854 


12 


2717 


-4.804 


-4.784 


0.1396 


0.0320 


0.5847 


0.000 


37.2 ± 5.4 


4.59 ± 0.53 


HD36108 


13 


2971 


-5.011 


-5.013 


0.0148 


0.0114 


0.2170 


0.000 


19.3 ± 2.7 


3.54 ± 0.41 


HD36003 


61 


1495 


-4.894 


-4.901 


0.0984 


0.0703 


0.1893 


0.000 


45.7 ± 6.6 


5.54 ± 0.64 


HD36379 


46 


2342 


-4.994 


-4.994 


0.0153 


0.0134 


0.0267 


1.000 


15.3 ±2.5 


2.92 ± 0.34 


HIP26542 


6 


1135 


-4.833 


-4.817 


0.1217 


0.0481 


0.6010 


0.000 


40.5 ±6.1 


4.80 ± 0.55 


HD37986 


16 


2156 


-5.050 


-5.047 


0.0225 


0.0133 


0.0677 


0.750 


46.1 ±4.3 


8.43 ± 0.97 


HD38277 


6 


2714 


-5.024 


-5.025 


0.0147 


0.0056 


0.6010 


0.000 


26.2 ±3.1 


4.77 ± 0.55 


HD38973 


24 


2354 


-4.966 


-4.969 


0.0278 


0.0245 


0.0732 


0.333 


18.2 ±2.7 


3.18 ±0.37 


HD40105 


19 


1794 


-5.177 


-5.176 


0.0198 


0.0150 


0.0716 


0.500 


59.9 ± 5.0 


11.53 ± 1.33 


HD40397 


11 


2756 


-5.024 


-5.023 


0.0190 


0.0177 


0.3235 


0.000 


36.5 ± 3.7 


6.65 ± 0.77 
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Name 


N 

1 Y meas 


Span 


10g^ K 


log^ K 


^hk> 


^HK> 


"■Urn 


Prob 


"rot 


Age 






[days] 


mean 


median 


raw 


short-term 




(A < 0.04) 


[days] 


[Gyr] 


HD44447 


19 


2745 


-4.995 


-4.992 


0.0319 


0.0252 


0.1234 


0.333 


15.3 ±2.5 


2.93 ± 0.34 


HD44120 


11 


2745 


-5.050 


-5.052 


0.0174 


0.0119 


0.1082 


0.250 


20.3 ± 2.7 


3.87 ± 0.45 


HD44594 


18 


2360 


-4.991 


-4.995 


0.0424 


0.0212 


0.1359 


0.417 


27.7 ± 3.2 


4.88 ± 0.56 


HD44573 


26 


2354 


-4.597 


-4.586 


0.1664 


0.0942 


0.5963 


0.000 


23.0 ± 5.2 


2.04 ± 0.23 


HD44420 


14 


2356 


-4.999 


-4.990 


0.0443 


0.0196 


0.2541 


0.000 


31.8 ±3.5 


5.64 ± 0.65 


HD45289 


12 


2749 


-5.033 


-5.034 


0.0201 


0.0190 


0.0936 


0.250 


30.6 ± 3.3 


5.62 ± 0.65 


HD47186 


98 


2635 


-5.023 


-5.021 


0.0146 


0.0132 


0.0143 


1.000 


35.3 ±3.6 


6.42 ± 0.74 


HD48611 


11 


1189 


-4.817 


-4.821 


0.0984 


0.0416 


0.5928 


0.000 


28.5 ±3.8 


4.08 ± 0.47 


HD50590 


15 


2333 


-4.999 


-5.001 


0.0797 


0.0343 


0.3572 


0.250 


50.9 ± 5.6 


7.68 ± 0.88 


HD50806 


45 


2084 


-5.095 


-5.093 


0.0199 


0.0169 


0.0318 


1.000 


38.6 ±3.6 


7.51 ±0.87 


HD52919 


7 


1067 


-4.772 


-4.770 


0.0346 


0.0232 


0.6010 


0.000 


36.7 ± 6.4 


3.87 ± 0.45 


HD59711A 


15 


2691 


-4.959 


-4.961 


0.0301 


0.0221 


0.1274 


0.250 


24.1 ±3.1 


4.09 ± 0.47 


HD65562 


18 


2332 


-4.954 


-4.953 


0.0884 


0.0486 


0.2262 


0.000 


44.7 ± 4.8 


7.08 ± 0.82 


HD66221 


19 


2351 


-5.034 


-5.036 


0.0263 


0.0199 


0.1183 


0.250 


37.0 ±3.7 


6.81 ± 0.78 


HD69655 


18 


2753 


-4.959 


-4.959 


0.0175 


0.0155 


0.0447 


0.917 


16.9 ± 2.6 


2.97 ± 0.34 


HD70889 


17 


1580 


-4.781 


-4.788 


0.1065 


0.0622 


0.3947 


0.333 


14.4 ± 2.8 


1.98 ±0.23 


HD71334 


16 


2355 


-4.995 


-4.993 


0.0188 


0.0164 


0.0935 


0.000 


25.3 ±3.1 


4.46 ±0.51 


HD72579 


18 


2707 


-5.062 


-5.062 


0.0229 


0.0181 


0.0734 


0.583 


45.9 ± 4.2 


8.52 ± 0.98 


HD72769 


22 


2704 


-5.052 


-5.052 


0.0210 


0.0143 


0.0574 


0.500 


40.2 ± 3.8 


7.51 ±0.86 


HD74014 


15 


2702 


-5.043 


-5.042 


0.0179 


0.0127 


0.1276 


0.333 


42.0 ± 4.0 


7.73 ± 0.89 


HD78429 


59 


1961 


-4.921 


-4.926 


0.0760 


0.0703 


0.1244 


0.000 


25.1 ±3.2 


4.09 ± 0.47 


HD78558 


15 


2354 


-5.020 


-5.020 


0.0158 


0.0109 


0.0499 


0.833 


23.5 ± 2.9 


4.26 ± 0.49 


HD78612 


19 


2703 


-5.030 


-5.030 


0.0256 


0.0232 


0.0733 


0.583 


22.4 ± 2.9 


4.13 ±0.48 


HD81639 


19 


2707 


-5.007 


-5.006 


0.0284 


0.0195 


0.1272 


0.167 


35.8 ±3.7 


6.40 ± 0.74 


HD83529 


23 


2707 


-4.993 


-4.993 


0.0148 


0.0122 


0.0405 


0.917 


17.6 ±2.6 


3.21 ± 0.37 


HD86171 


8 


1209 


-4.807 


-4.807 


0.1104 


0.0328 


0.6010 


0.000 


28.8 ±3.9 


4.05 ± 0.47 


HD88084 


18 


2693 


-4.982 


-4.981 


0.0157 


0.0099 


0.0918 


0.833 


26.1 ±3.2 


4.55 ± 0.52 


HD88218 


19 


2701 


-5.055 


-5.056 


0.0123 


0.0089 


0.0909 


0.833 


23.1 ±2.9 


4.36 ± 0.50 


HD88742 


23 


1869 


-4.694 


-4.696 


0.2041 


0.1354 


0.5632 


0.000 


11.3 ±2.7 


1.37 ±0.16 


HD92588 


33 


1883 


-5.167 


-5.167 


0.0216 


0.0185 


0.0424 


0.917 


58.7 ±4.9 


11.30 ± 1.30 


HD92719 


15 


2324 


-4.860 


-4.861 


0.0668 


0.0500 


0.4446 


0.000 


18.6 ±2.9 


2.83 ± 0.33 


HD93385 


121 


2694 


-4.984 


-4.984 


0.0147 


0.0135 


0.0187 


1.000 


18.6 ±2.7 


3.32 ± 0.38 


HD95456 


77 


2339 


-4.963 


-4.965 


0.0447 


0.0318 


0.0607 


0.417 


11.6 ±2.3 


2.54 ± 0.29 


HD95521 


15 


2349 


-4.886 


-4.905 


0.1364 


0.0497 


0.5493 


0.000 


21.6 ±3.1 


3.38 ± 0.39 


HD96423 


21 


2703 


-5.021 


-5.023 


0.0189 


0.0158 


0.0619 


0.750 


31.4 ±3.4 


5.70 ± 0.66 


HD97037 


18 


2669 


-5.004 


-5.003 


0.0139 


0.0117 


0.1007 


0.583 


21.7 ±2.9 


3.89 ± 0.45 


HD97343 


20 


2677 


-5.021 


-5.022 


0.0169 


0.0133 


0.0605 


0.750 


40.8 ± 4.0 


7.36 ± 0.85 


HD97998 


17 


2352 


-4.948 


-4.948 


0.0164 


0.0144 


0.0958 


0.583 


22.5 ± 3.0 


3.79 ± 0.44 


HD98356 


8 


2326 


-4.786 


-4.790 


0.0905 


0.0445 


0.6010 


0.000 


32.5 ± 4.5 


4.30 ± 0.49 


HD102117 


32 


2282 


-5.053 


-5.053 


0.0240 


0.0224 


0.0440 


0.833 


37.9 ±3.7 


7.10 ±0.82 


HD 10243 8 


39 


2244 


-4.982 


-4.982 


0.0134 


0.0127 


0.0211 


1.000 


29.9 ± 3.4 


5.21 ± 0.60 


HD 104006 


26 


2233 


-5.054 


-5.053 


0.0208 


0.0177 


0.0572 


0.583 


44.5 ±4.1 


8.24 ± 0.95 


HD 104067 


83 


2271 


-4.752 


-4.753 


0.0975 


0.0872 


0.1218 


0.000 


34.0 ± 5.6 


3.82 ± 0.44 


HD 104263 


29 


2245 


-5.035 


-5.033 


0.0284 


0.0263 


0.0725 


0.500 


40.5 ± 3.9 


7.42 ± 0.85 


HD 104982 


39 


2271 


-4.974 


-4.973 


0.0140 


0.0133 


0.0316 


1.000 


25.8 ± 3.2 


4.46 ±0.51 


HD105837 


21 


2652 


-4.877 


-4.879 


0.0558 


0.0360 


0.1903 


0.417 


13.9 ±2.6 


2.27 ± 0.26 


HD 109409 


22 


2649 


-5.110 


-5.111 


0.0207 


0.0145 


0.0570 


0.833 


35.1 ±3.4 


6.92 ± 0.80 


HD1 10619 


17 


2648 


-4.922 


-4.925 


0.0527 


0.0254 


0.2233 


0.250 


25.1 ±3.2 


4.10 ±0.47 


HD111031 


26 


2245 


-5.036 


-5.036 


0.0176 


0.0146 


0.0479 


0.833 


33.4 ± 3.5 


6.14 ±0.71 


HD 112540 


7 


2295 


-4.774 


-4.771 


0.2205 


0.0781 


0.6010 


0.000 


24.8 ± 3.7 


3.35 ± 0.39 


HD1 14853 


37 


3011 


-4.957 


-4.959 


0.0393 


0.0376 


0.0842 


0.083 


24.0 ±3.1 


4.07 ± 0.47 


HD1 15585 


19 


2641 


-5.076 


-5.084 


0.0299 


0.0248 


0.1113 


0.250 


41.4 ±3.8 


7.90 ± 0.91 


HD 116920 


8 


2305 


-4.827 


-4.820 


0.1210 


0.0834 


0.6010 


0.000 


37.9 ±5.1 


4.99 ± 0.57 


HD117105 


18 


2658 


-4.976 


-4.976 


0.0128 


0.0113 


0.0529 


0.750 


17.3 ±2.6 


3.09 ± 0.36 


HD 117207 


15 


2642 


-5.035 


-5.035 


0.0076 


0.0069 


0.0985 


0.667 


37.0 ±3.7 


6.82 ± 0.78 


HD1 19782 


17 


2636 


-4.697 


-4.703 


0.0945 


0.0811 


0.4908 


0.000 


28.2 ± 4.7 


3.18 ±0.37 


HD123265 


17 


2642 


-5.073 


-5.078 


0.0250 


0.0207 


0.1342 


0.000 


49.8 ± 4.5 


9.14 ± 1.05 


HD122862 


18 


2295 


-5.027 


-5.026 


0.0129 


0.0113 


0.0647 


0.667 


18.5 ±2.6 


3.48 ± 0.40 


HD 124292 


25 


3006 


-5.009 


-5.012 


0.0243 


0.0208 


0.0791 


0.083 


37.0 ±3.8 


6.62 ± 0.76 


HD124364 


14 


2636 


-4.865 


-4.880 


0.1139 


0.0357 


0.5563 


0.000 


24.1 ±3.3 


3.68 ± 0.42 


HD125184 


68 


2131 


-5.029 


-5.025 


0.0411 


0.0346 


0.0880 


0.500 


36.8 ±3.7 


6.72 ± 0.77 


HD125455 


18 


2165 


-4.916 


-4.915 


0.1211 


0.0543 


0.5083 


0.250 


42.3 ± 4.8 


6.43 ± 0.74 


HD126525 


42 


2518 


-4.993 


-4.993 


0.0138 


0.0118 


0.0308 


1.000 


30.3 ± 3.4 


5.34 ± 0.62 


HD 128674 


19 


2263 


-4.958 


-4.957 


0.0188 


0.0129 


0.0809 


0.500 


27.7 ± 3.3 


4.71 ±0.54 


HD 130992 


18 


1831 


-4.858 


-4.860 


0.0722 


0.0505 


0.3187 


0.000 


42.2 ±6.1 


5.17 ±0.59 


HD130930 


18 


2117 


-5.016 


-5.015 


0.0407 


0.0241 


0.1726 


0.000 


50.3 ±5.1 


8.22 ± 0.95 
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Name 


N 

1 Y meas 


Span 


10g^ K 


log^ K 


^hk> 


^hk) 


"■Urn 


Prob 


"rot 


Age 






[days] 


mean 


median 


raw 


short-term 




(A < 0.04) 


[days] 


[Gyr] 


HD 134060 


100 


2616 


-4.983 


-4.983 


0.0172 


0.0158 


0.0209 


1.000 


22.3 ± 2.9 


3.91 ± 0.45 


HD134985 


15 


2619 


-5.068 


-5.074 


0.0309 


0.0217 


0.3366 


0.000 


44.3 ±4.1 


8.34 ± 0.96 


HD 134606 


107 


2489 


-5.049 


-5.049 


0.0138 


0.0133 


0.0155 


1.000 


40.1 ±3.8 


7.46 ± 0.86 


HD136894 


38 


2044 


-5.006 


-5.007 


0.0130 


0.0118 


0.0267 


1.000 


35.7 ±3.7 


6.39 ± 0.74 


HD137303 


16 


2044 


-4.811 


-4.807 


0.1328 


0.0574 


0.4776 


0.000 


39.0 ±6.1 


4.48 ± 0.52 


HD140901 


24 


1546 


-4.694 


-4.697 


0.1380 


0.0879 


0.5881 


0.000 


20.3 ± 3.6 


2.43 ± 0.28 


HD 142709 


13 


2117 


-5.049 


-5.058 


0.1120 


0.0406 


0.4298 


0.000 


56.9 ± 6.7 


8.08 ± 0.93 


HD143114 


19 


1789 


-4.991 


-4.989 


0.0115 


0.0094 


0.0734 


0.667 


21.3 ±2.9 


3.77 ± 0.43 


HD 1444 11 


11 


1719 


-4.851 


-4.843 


0.0709 


0.0541 


0.5024 


0.000 


40.9 ± 5.6 


5.22 ± 0.60 


HD 144585 


14 


2539 


-5.033 


-5.043 


0.0449 


0.0098 


0.4071 


0.000 


29.3 ± 3.2 


5.37 ± 0.62 


HD145809 


13 


2174 


-5.062 


-5.061 


0.0115 


0.0091 


0.0704 


0.750 


24.8 ± 2.9 


4.68 ± 0.54 


HD145598 


36 


2107 


-5.004 


-5.003 


0.0256 


0.0222 


0.0584 


0.833 


28.2 ± 3.2 


5.03 ± 0.58 


HD 145666 


22 


1410 


-4.778 


-4.777 


0.0543 


0.0390 


0.2174 


0.083 


14.3 ± 2.8 


1.96 ± 0.23 


HD147512 


29 


1734 


-4.999 


-4.999 


0.0126 


0.0103 


0.0322 


1.000 


35.4 ± 3.7 


6.28 ± 0.72 


HD150433 


50 


2077 


-5.000 


-5.000 


0.0107 


0.0097 


0.0181 


1.000 


24.2 ± 3.0 


4.29 ± 0.49 


HD151504 


14 


1494 


-5.031 


-5.032 


0.0097 


0.0081 


0.0629 


0.667 


42.3 ±4.1 


7.68 ± 0.88 


HD154363 


17 


2516 


-4.906 


-4.910 


0.1488 


0.0723 


0.4597 


0.000 


46.4 ± 6.5 


5.75 ± 0.66 


HD157347 


20 


2893 


-5.011 


-5.010 


0.0144 


0.0095 


0.0684 


0.667 


31.0 ±3.4 


5.57 ± 0.64 


HD157338 


25 


1461 


-4.975 


-4.972 


0.0232 


0.0189 


0.0850 


0.333 


18.4 ±2.7 


3.25 ± 0.37 


HD161612 


32 


2154 


-5.014 


-5.014 


0.0117 


0.0107 


0.0221 


1.000 


39.4 ± 3.9 


7.07 ± 0.81 


HD162236 


16 


1069 


-4.709 


-4.714 


0.1443 


0.1302 


0.6010 


0.000 


22.4 ± 3.8 


2.74 ± 0.32 


HD 162396 


39 


1885 


-5.007 


-5.006 


0.0203 


0.0170 


0.0414 


0.917 


11.3 ±2.2 


2.94 ± 0.34 


HD165920 


13 


1791 


-5.048 


-5.048 


0.0239 


0.0185 


0.1102 


0.000 


48.9 ± 4.6 


8.72 ± 1.00 


HD 166724 


47 


2526 


-4.740 


-4.739 


0.1268 


0.0958 


0.2286 


0.000 


30.9 ± 4.7 


3.74 ± 0.43 


HD168871 


22 


2537 


-4.988 


-4.989 


0.0180 


0.0145 


0.0500 


0.833 


18.8 ±2.7 


3.36 ± 0.39 


HD 170493 


13 


2101 


-4.802 


-4.802 


0.2087 


0.0873 


0.6010 


0.000 


38.7 ±6.3 


4.29 ± 0.49 


HD171665 


13 


1730 


-4.914 


-4.906 


0.0505 


0.0149 


0.2931 


0.000 


28.2 ± 3.5 


4.57 ± 0.53 


HD171990 


34 


1797 


-5.101 


-5.099 


0.0237 


0.0213 


0.0594 


0.667 


21.6 ±2.7 


4.31 ±0.50 


HD 174545 


13 


2535 


-4.911 


-4.894 


0.1092 


0.0521 


0.5910 


0.000 


42.4 ± 4.9 


6.35 ± 0.73 


HD176986 


80 


2474 


-4.833 


-4.829 


0.0790 


0.0747 


0.1288 


0.250 


39.0 ± 5.3 


5.04 ± 0.58 


HD177758 


7 


2749 


-4.991 


-4.993 


0.0087 


0.0054 


0.5928 


0.000 


16.4 ± 2.6 


3.04 ± 0.35 


HD 177409 


18 


1980 


-4.867 


-4.878 


0.0866 


0.0430 


0.2814 


0.083 


16.7 ±2.8 


2.58 ± 0.30 


HD 180409 


8 


2747 


-4.942 


-4.942 


0.0094 


0.0044 


0.3603 


0.000 


15.3 ±2.6 


2.70 ±0.31 


HD181433 


154 


2882 


-5.100 


-5.100 


0.0276 


0.0261 


0.0211 


1.000 


58.4 ±5.8 


9.47 ± 1.09 


HD183658 


14 


1695 


-4.983 


-4.980 


0.0185 


0.0120 


0.0875 


0.750 


24.9 ±3.1 


4.34 ± 0.50 


HD183783 


8 


2526 


-4.954 


-4.942 


0.1540 


0.0418 


0.6010 


0.000 


49.7 ± 6.5 


6.54 ± 0.75 


HD185615 


22 


2305 


-5.033 


-5.032 


0.0294 


0.0259 


0.0795 


0.000 


38.1 ±3.8 


7.00 ± 0.81 


HD 187456 


11 


2145 


-4.925 


-4.926 


0.0304 


0.0144 


0.2144 


0.417 


47.1 ±6.2 


6.14 ±0.71 


HD 189625 


20 


1751 


-4.788 


-4.779 


0.1556 


0.0742 


0.4726 


0.250 


19.2 ±3.2 


2.65 ± 0.30 


HD 188748 


19 


2122 


-4.981 


-4.980 


0.0209 


0.0183 


0.0732 


0.333 


28.6 ± 3.3 


4.98 ± 0.57 


HD 189567 


154 


2742 


-4.941 


-4.943 


0.0390 


0.0344 


0.0421 


0.917 


24.7 ± 3.2 


4.11 ±0.47 


HD 190248 


84 


2491 


-5.052 


-5.052 


0.0164 


0.0137 


0.0244 


1.000 


41.4 ±3.9 


7.70 ± 0.89 


HD 190954 


9 


2111 


-5.015 


-5.017 


0.0279 


0.0116 


0.4987 


0.000 


35.0 ±3.6 


6.31 ±0.73 


HD 192031 


8 


2525 


-5.051 


-5.051 


0.0107 


0.0067 


0.4282 


0.000 


39.0 ±3.8 


7.29 ± 0.84 


HD193193 


31 


1760 


-4.938 


-4.943 


0.0411 


0.0313 


0.1614 


0.000 


17.5 ±2.7 


2.96 ± 0.34 


HD195564 


11 


2731 


-5.106 


-5.104 


0.0231 


0.0098 


0.3661 


0.000 


36.3 ± 3.5 


7.15 ±0.82 


HD196761 


8 


2730 


-4.951 


-4.949 


0.0667 


0.0180 


0.5672 


0.000 


33.1 ±3.7 


5.58 ± 0.64 


HD197210 


7 


2073 


-4.890 


-4.880 


0.0933 


0.0098 


0.5845 


0.000 


29.3 ± 3.6 


4.61 ± 0.53 


HD199288 


14 


2491 


-4.963 


-4.962 


0.0101 


0.0052 


0.0658 


0.667 


18.1 ±2.7 


3.16 ±0.36 


HD203384 


9 


2333 


-5.002 


-5.001 


0.0489 


0.0175 


0.5404 


0.000 


39.8 ±4.0 


7.04 ±0.81 


HD203850 


7 


2871 


-4.921 


-4.924 


0.0534 


0.0159 


0.6010 


0.000 


44.3 ± 5.2 


6.49 ± 0.75 


HD204385 


7 


2222 


-4.967 


-4.968 


0.0280 


0.0148 


0.6010 


0.000 


19.4 ± 2.8 


3.37 ± 0.39 


HD205536 


23 


2219 


-5.026 


-5.022 


0.0397 


0.0337 


0.2436 


0.000 


41.1 ±4.0 


7.44 ± 0.86 


HD207700 


15 


2111 


-4.999 


-5.000 


0.0245 


0.0191 


0.0897 


0.000 


33.0 ±3.5 


5.87 ± 0.68 


HD209742 


7 


1069 


-4.872 


-4.891 


0.1331 


0.0054 


0.6005 


0.000 


39.1 ±4.7 


5.68 ± 0.65 


HD210752 


12 


1779 


-4.934 


-4.934 


0.0134 


0.0091 


0.1771 


0.167 


12.1 ±2.3 


2.36 ± 0.27 


HD211038 


73 


1880 


-5.214 


-5.214 


0.0086 


0.0079 


0.0153 


1.000 


62.3 ± 5.0 


12.35 ± 1.42 


HD210918 


31 


1858 


-5.010 


-5.008 


0.0249 


0.0244 


0.0304 


1.000 


27.1 ±3.2 


4.86 ± 0.56 


HD212708 


34 


1072 


-5.046 


-5.047 


0.0233 


0.0201 


0.0731 


0.583 


38.8 ±3.8 


7.21 ±0.83 


HD2 13628 


7 


2176 


-4.953 


-4.951 


0.0328 


0.0052 


0.6010 


0.000 


33.2 ±3.7 


5.62 ± 0.65 


HD213941 


26 


2218 


-4.961 


-4.961 


0.0333 


0.0304 


0.1169 


0.333 


27.8 ± 3.3 


4.74 ± 0.55 


HD215456 


94 


2616 


-5.089 


-5.087 


0.0206 


0.0185 


0.0224 


1.000 


28.5 ±3.1 


5.51 ±0.63 


HD2 19249 


22 


1283 


-4.960 


-4.959 


0.0283 


0.0168 


0.1333 


0.000 


30.2 ± 3.5 


5.15 ±0.59 


HD220256 


23 


2222 


-5.033 


-5.034 


0.0320 


0.0222 


0.1114 


0.500 


48.6 ± 4.7 


8.49 ± 0.98 


HD220339 


11 


2122 


-4.836 


-4.818 


0.1627 


0.0379 


0.5919 


0.000 


37.6 ±4.9 


5.13 ±0.59 


HD220507 


50 


2221 


-5.050 


-5.051 


0.0182 


0.0149 


0.0455 


0.917 


34.0 ± 3.5 


6.34 ± 0.73 
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Table 3. Continued. 



Name 


N 

1 ¥ meas 


Span 


10g^ K 


log^ K 


^hk> 


^HK> 


"■Urn 


Prob 


"rot 


Age 






[days] 


mean 


median 


raw 


short-term 




(A < 0.04) 


[days] 


[Gyr] 


HD221356 


15 


2165 


-4.940 


-4.939 


0.0329 


0.0251 


0.232A 


0.000 


11.2 ±2.3 


2.40 ± 0.28 


HD221420 


28 


1420 


-5.158 


-5.155 


0.0244 


0.0204 


0.0641 


0.667 


37.0 ± 3.4 


7.62 ± 0.88 


HD222422 


7 


2113 


-4.792 


-4.791 


0.0671 


0.0365 


0.6010 


0.000 


26.4 ± 3.8 


3.66 ± 0.42 


HD223121 


13 


2133 


-4.845 


-4.848 


0.1415 


0.0867 


0.6010 


0.000 


39.8 ±5.3 


5.22 ± 0.60 


HD223171 


64 


1809 


-4.995 


-4.995 


0.0476 


0.0439 


0.0695 


0.333 


27.9 ± 3.2 


4.92 ± 0.57 


HD224619 


11 


2136 


-5.000 


-5.001 


0.0256 


0.0091 


0.1550 


0.417 


37.6 ±3.8 


6.67 ± 0.77 
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